UNCLASSIFIED 


w 


kepAoduced 
luf  the 


ARMED  SERVICES  TECHNICAL  INFORMAITON  AGENCY 
ARLINGTON  HALL  STATION 
ARLINCIW  12.  VIRGINIA 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S» 
Government  thereby  Incurs  no  responsibility,  nor  any 
obligation  lAatsoeverj  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  maimer  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rlj^ts 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


swc 

TR 

61-67 


SWC-TR-61-67 


iO 

HEADQUARTERS 

AIR  FORCE  SPECIAL  WEAPONS  CENTER 
AIR  FORCE  SYSTEMS  COMMAND 


A  THEORETICAL  STUDY  OF  RESPONSE  OF  SOLIDS 
TO  IMPULSIVE  LOADS  OF  HIGH  PRESSURE 

by 

Joseph  J.  Poczatek 
L.  E.  F'ugelso 

Final  Report 
-October  1961 


American  Machine  &  Foundry  Company 
Mechanics  Research  Division 
Niles,  Illinois 


This  research  is  a  part 
of  Project  DEFENDER 
sponsored  by  the 
Advanced  Research 
Projects  Agency 
Department  of  Defense 


ftSTIA 


HEADQUARTERS 

AIR  FORCE  SPECIAL  WEAPONS  CENTER 
Air  Force  Systems  Command 
Kirtland  Air  Force  Base 
New  Mexico 


When  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  in  connection  with  a  definitely  related  Government 
procurement  operation,  the  United  States  Government  thereby  Incurs  no 
responsibility  nor  any  obligation  whatsoever;  and  the  fact  that  the  Government 
may  have  formulated,  furnished,  or  in  any  way  supplied  the  said  drawings, 
specifications,  or  other  data,  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any  other  person  or  corporation, 
or  conveying  any  rights  or  permission  to  manufacture,  use,  or  sell  any 
patented  invention  that  may  in  any  way  be  related  thereto. 

This  report  is  made  available  for  study  upon  the  understanding  that  the 
Government's  proprietary  interests  in  and  relating  thereto  shall  not  be  im¬ 
paired,  In  case  of  apparent  conflict  between  the  Government's  proprietary 
interests  and  those  of  others,  notify  the  Staff  Judge  Advocate,  Air  Force 
Systems  Command,  Andi-ews  AF  Base,  Washington  Z5,  DC. 


This  report  is  published  for  the  exchange  and  stimulation  of  Ideas;  it  does 
not  necessarily  express  the  intent  or  policy  of  any  higher  headquarters. 

Qualified  requesters  may  obtain  copies  of  this  report  from  ASTIA. 

Orders  will  be  expedited  if  placed  through  the  librarian  or  other  staff  member 
designated  to  request  and  receive  documents  from  ASTIA, 


I 

I 


AFSWC-TR-61-67 


A  THEORETICAL  STUDY  OF  RESPONSE  OF  SOLIDS 
TO  IMPULSIVE  LOADS  OF  HIGH  PRESSURE 


by 

Joseph  J.  Poczatek 
L.  E.  Fugelso 

Final  Report 
October  1961 


American  Machine  &  Foundry  Company 
Mechanics  Research  Division 
Niles,  Illinois 


Research  Directorate 
AIR  FORCE  SPECIAL  WEAPONS  CENTER 
Air  Force  Systems  Command 
Kirtland  Air  Force  Base 
New  Mexico 


Approved: 


Director,  Research  Directorate 


Project  4778 

ARPA  Order  6-5  8,  Task  2Z 
Contract  AF  29(601)-2533 
AMF-MRD  Project  MR  1134 


I 


FOREWORD 

The  Mechanics  Research  Division  of  American  Machine  and  Foundry  Con5)any 
takes  pleasure  in  presenting  this  technical  report  to  the  Physics  D.!  v’i.  2 i  .;n , 

Re  Gear  oh  Direct. .■)  rate,,  of  the  Air  Force  Special  Weapons  Center^  Kii'tland  AFB, 

This  technical  report  contains  the  results  of  theoretical  research  on 
plastic  deformation  vmder  high  short  duration  loading.  The  research  was  carried 
out  as  a  portion  of  Contract  AF  29(601 )-2533. 

The  remaining  effort  on  this  contract  has  resulted  in  a  separate  report, 
designated  TIR-62-5,  which  treats  the  deformational  behavior  of  solids  under 
high  iii5)ulsive- pressures  of  significantly  shorter  time  durations. 

Acknowledgement  is  made  to  Mr.  Joseph  J.  Poczatek,  the  project  supervisor, 
and  Dr.  G.L.  Neidhardt,  both  of  American  Machine  and  Foundry  Con^jany,  Mechanics 
^search  Dlvlsioii  and  Dr,  S.  Raynor,  Northwestern  University,  Consultant  to 
Mechanics  Research  Division,  for  stimulating  conversations  and  discussions, 

J.D.  Stein,  H-M.  Gruen,  M.F.  Darienzo  and  D.H.  Duel  i>erformed  the  numerical 
computation,  drew  the  graphs  and  proof read  the  report .  Capt.  Marvin  Atkins  of' 
the  Physic.G  Division,  AFSV7C,  was  the  technical  raoniter.  The  principal  inves¬ 
tigator  on  this  jxjrtlon  of  the  project  was  L.E.  Fxjgelso. 


Respectfully  submitted, 

APPROVED:  AMERICAN  MACHINE  &  FOUNDRY  C(»lPAHr 

Mechanics  Research  Division 


G.L.  Neidhardt,  Group  Leader 
Engineering  Mathematics  Group 
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ABSTRAC  T 


The  response  of  materials  and  structures  to  pressures  on  the  order  of  one 
megabar  is  considered.  The  problem  is  decomposed  according  to  the  charac¬ 
teristic  response  time  for  various  failure  modes.  The  short-time  effects 
study  analyzes  propagation  of  these  high-pressure  pulses.  General  equations 
are  developed  for  a  solid  under  finite  strain,  with  dissipation  mechanism, 
heat  conduction,  and  elasticity  coefficients  open  to  specification.  Integral 
equations  following  motion  of  non-isentropic  waves  are  obtained  for  uni-axial 
strain.  Conditions  for  instantaneously  forming  finite  discontinuities  are  de¬ 
lineated.  Succeeding  states  are  sought  isentropically,  resulting  in  a  theo¬ 
retical  isentrope.  A  closed  form  family  of  isentropes  is  obtained  for  aluminum 
using  the  Hugoniot. 

The  long-time  effects  study  develops  dynamic  elastm-plastic  deformation 
theory  of  solids  utilizing  dislocation  theory.  Continuum  equations  for  the 
elastic-plastic  displacements  and  stresses  are  derived.  Propagation  of  the 
plastic  wave,  apparent  time-dependent  increase  in  yield  conditions,  and  other 
deformational  characteristics  are  discussed.  The  elastic-plastic  deformation 
equations  are  solved  for  three  particular  cases; 

1.  Propagation  of  a  stress  wave  in  a  one-dimensional  rod. 

2.  Deformation  of  a  semi-infinite  half  space  under  a  suddenly  applied 
strip  load. 

3.  Deformation  of  an  infinitely  long  thin  cylindrical  shell  under  a  suddenly 
applied  strip  load  along  a  generator. 
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PLASTIC  DEFORMATION  UNDER  IMPACT  LOADING 


A  theoretical  approach  to  the  dynamic  elastic-plastic  deformation  of 
poly ciy stall ine  solids  is  developed  utilizing  the  theories  of  dislocations  and 
irreversible  thermodynamics. 

In  a  solid  body,  stresses  which  are  above  the  static  yield  condition  may  be 
developed  under  time  dependent  loading.  If  a  state  of  stress  which  falls  outside 
the  yield  conditions  comes  into  existence,  the  solid  body  will  react  to  alleviate 
this  imstable  state  by  the  mechanism  of  moving  dislocations. 

Equations  for  the  dynamic  elastic-plastic  displacements  and  stresses  are 
derived  on  the  basis  of  a  model  for  a  continuum  density  of  dislocations.  This 
model  for  elastic-plastic  deformation  is  established  by  considering  the  dislocations 
as  being  randomly  distributed.  The  physical  properties  of  the  dislocations  and 
their  reaction  to  applied  stresses  are  derived  following  the  basic  geometrical 
theories  of  Kondo,  Kroner,  Taylor,  and  Bilby. 

The  elastic-plastic  deformation  equations  are  solved  for  three  particular 
cases  of  interest.  The  first  is  the  one- dimensional  rod  impacted  at  one  end. 

The  second  is  the  elastic-plastic  deformation  of  a  semi-infinite  half- space  under 
a  suddenly  applied  strip  load.  The  dynamic  plastic  deformation  of  a  cylindrical 
shell  is  the  third  problem  studied.  Numerical  solutions  for  the  half-space 
and  the  shell  problems  are  given. 

Several  featiires  of  elastic-plastic  stress  wave  propagation  and  stress, 

states  under  time  dependent  loads  are  discussed  in  some  detail.  In  particiilar, 

the  propagation  of  the  so-called  plastic  wave,  the  apparent  increase  in  yield 
conditions  under  time  dependent  loads,  and  the  course  of  deformation  are  discussed. 
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I.  Introduction 


When  a  material  specimen  is  subjected  to  external  loads,  it  deforms. 

If  the  loads  are  sufficiently  small  the  resulting  deformation  is  elastic, 
i.e.,  the  deformation  gradients  are  almost  linear  functions  of  the  stresses 
developed  Internally.  When  the  loads  are  higher  than  a  specified  level,  the 
defonnation  ceases  to  be  perfectly  elastic.  Normally  the  deformation  increases 
more  rapidly  than  the  stress.  The  basic  physical  feature  of  the  deformation 
is  that  the  material  slips  along  certain  lines  or  planes  within  the  body. 

This  deformation  is  plastic  deformation  of  the  medium. 

The  problem  of  plastic  deformation  arises  in  many  areas  of  engineering 
and  applied  science.  Various  in^^et  problems,  where  the  peak  stresses  rise 
far  above  the  proportional  limit,  occur  repeatedly;  for  example,  the  deforma¬ 
tion  of  a  plate  under  iuq^ct  from  a  bullet.  In  seismology  and  the  study  of 
earthquakes,  the  deformation  close  to  the  focus  of  the  earthquake  is  certainly 
in  the  plastic  region. 

/Q\  1  / 

The  classical  approach  to  the  static  problem'  '  —  is  that  there  exists  some 
relation  among  the  stresses  which  cannot  be  exceeded.  At  incipient  plastic 
deformation,  the  state  of  stress  is  at  the  maximum  point  or  locus  of  this 
relation,  called  the  yield  condition.  Insertion  of  this  condition  into  the 
equation  of  equilibrliun  of  an  elastic  solid  yields  a  solution  for  the  stress 
Internal  to  the  body. 

When  the  stress  field  is  caused  by  a  load  applied  as  a  function  of  time, 
many  difficulties  ensue.  This  problem  of  dynamic  elastic-plastic  deformation 
has  been  treated  by  several  authors, (l0),(ll),(2l), (22)(23)  approach 

Superscripts  refer  to  references  listed  in  the  bibliography. 


1 


to  the  dynamic  plastic  deformation  is  to  assume  a  non-linear  stress-strain^^^^' 

(ll)  (22) 

'  ^  curve  and  insert  this  expression  into  the  equations  of  motion  for 

an  ela.stic  body.  Several  difficulties  arise  immediately.  First  the  equations 

are  non-linear  and  do  not  lend  themselves  to  easy  solution.  Second,  as  pointed 

(21) 

out  by  Taylor,  the  non-linear  stress-strain  curve  is  a  function  of  the  rate 
of  loading. 

In  either  approach  to  the  problem  of  inelastic  deformation,  certain  in¬ 
adequacies  are  readily  apparent.  Particularly  when  the  loads  are  applied 
over  a  very  short  time  duration  or  applied  very  rapidly,  the  experimental 
evidence  does  not  agree  closely  with  the  theoretical  results.  Under  rapidly 
applied  loads,  the  yield  point  is  raised  and  the  deformation  itself  becomes 
a  function  of  strain  rate  or  loading  rate,^^^’  (^5) 

To  overcome  these  difficulties,  the  deformation  of  a  solid  body  under 
stresses  which  are  in  the  region  above  the  elastic  stress  stability  region 
Is  studied  by  examining  the  mechanism  of  deformation. 

Above  the  yield  stress,  experimental  evidence  shows  that  for  many  metals, 
the  body  has  \mdergone  deformation  by  slipping  along  certain  lines  or  planes. 

This  slip  mechanism  for  plastic  deformation  has  been  postulated  as  the 
basic  physical  phenomenon  imderlying  plastic  deformation.  In  a  time  dependent 
slip,  it  is  unlikely  that  all  the  atoms  along  the  slip  plane  move  at  once. 

The  energy  requirements  for  this  are  many  times  the  actual  energy  involved. 

It  becomes  apparent  that  the  slip  occurs  a  little  bit  at  a  time  to  account 
for  the  energy  involved. 

In  193^  Taylor, Orowan, and  Polyani^^^^  introduced  the  concept 
of  plastic  deformation  by  the  movement  of  dislocations  and  iii5)erfections  in 
crystals.  This  basic  concept  remains  as  the  foundation  of  plastic  deformation. 
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These  papers  pointed  out  that  when  plastic  flow  occurs  it  is  very  un¬ 
likely  that  an  entire  plane  of  atoms  move  all  at  once,  hut  rather  the  slip 
proceeds  piecemeal  along  the  slip  plane.  At  the  line  in  the  slip  plane  where 
slip  has  occurred  on  one  side  hut  not  in  the  rest  of  the  slip  plane,  there  is 
an  in5)erfection  in  the  crystal  lattice. 

The  mechanism  of  plastic  deformation  may  he  explained  simply.  In  a  real 
crystal  there  exists  a  number  of  inperfections,  either  pre-existing  or  created 
during  deformation,  which  interrupt  the  perfect  crystal  lattice.  An  example 
of  this  type  of  imperfection  is  where  a  row  or  plane  of  atoms  is  effectively 
missing  from  the  perfect  lattice.  To  take  up  the  space  and  retain  the  crystal 
structure,  the  lattice  distorts  itself  to  accommodate  the  inperfectlon.  The 
mechanism  of  plastic  deformation  may  be  stated  qualitatively  then  as  the  forma¬ 
tion  of  these  types  of  inperfections  and/or  the  subsequent  movement  or  re¬ 
arrangement  of  the  dislocations. 

The  theory  of  plastic  deformation  by  dislocations  has  advanced  rapidly. 

Many  properties  of  plastic  deformation  and  other  phenomena  in  the  solid  state 
may  be  described  by  appeal  to  dislocations.  Among  these  are  diffusion  and 
precipitation  in  solids  and  some  surface  phenomena. 

A  qualitative  description  of  plastic  deformation  by  the  dislocation  move¬ 
ment  is  a  necessary  preliminary  to  further  understanding  and  use  of  dislocations. 
Consider  the  two-dimensional  picture  of  a  square  crystal  lattice.  Figure  1. 

Figure  la  represents  the  original  vuideformed  lattice.  Under  stress,  slip 
starts  at  some  point  in  the  lattice  and  proceeds,  one  atom  slipping  at  a  time, 
until  an  Intermediate  configuration  (Figure  lb)  is  reached.  The  deformation 
proceeds,  then,  atom  by  atom,  until  the  slip  is  complete  (Figure  Ic). 
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This  dislocation  is  called  an  edge  dislocation.  The  displacement  here 
is  perpendicular  to  the  line  of  the  dislocation.  Other  types  of  dislocations 
are  possible,  the  main  one  being  the  screw  dislocation,  in  which  the  displace¬ 
ment  in  the  ^lip  plane  is  parallel  to  the  edge  of  the  dislocation. 

Another  class  of  imperfections  which  may  be  conveniently  described  by 
dislocations  is  the  interface  between  two  crystals  in  a  polycrystalline  specimen. 
If  the  angle  between  two  different  crystals  is  small,  the  mismatch  at  the 
boundary  is  described  as  a  row  of  edge  dislocations.  Thus  in  a  polycrystal, 
a  natural  source  of  dislocations  is  present. 

In  the  theory  of  plastic  flow  developed  by  Taylor,  the  flow  is  due  to 
movement  and  formation  of  dislocations.  Treating  these  phases  one  at  a  time, 
the  movement  of  dislocations  will  first  be  considered. 

In  order  to  move  the  dislocation  a  force  must  be  applied  to  the  atom  on 
the  undisturbed  side  of  the  dislocation.  The  atom,  or  more  accurately  speaking, 
a  row  of  atoms  along  the  dislocation  line  must  be  forced  over  an  energy  barrier. 
For  a  perfect  crystal  the  stress  corresponding  to  this  is  the  stress  which  pro¬ 
duces  a  shear  strain  of  l/2.  If  a  perfect  crystal  with  but  a  single  vacancy  is 
considered,  the  strained  atom  must  still  have  a  strain  of  l/2  to  come  to  the 
position  of  neutral  equilibrium  between  the  two  possible  sites.  The  stress  that 
has  to  be  applied  to  the  atom  is  lower,  however,  as  the  vacancy  causes  a  distor¬ 
tion  of  the  lattice. 

Similarly,  a  dislocation  alters  the  strain  energy  aroiuid  the  dislocation. 

The  presence  of  a  dislocation  in  a  crystal  distorts  the  lattice  surrounding 

the  dislocation  and  alters  the  strain  energy  there.  A  quantitative  expression 

(5) 

for  the  stress  due  a  dislocation,  formulated  by  Cottrell  is'  ' 
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where  r,  q  are  polar  coordinates,  h  is  the  displacement  due  to  the  dislo¬ 
cation,  and  ju  is  the  elastic  shear  modulus,  f (  0)  is  the  angular  distribu¬ 
tion  and  depends  on  the  geometry  of  the  lattice  and  the  type  of  dislocation 
involved. 

The  explicit  form  of  the  stress  around  the  dislocation  is  an  involved 
problem  and  will  not  be  considered  here. 

The  stress  required  to  move  the  dislocation  then  is  of  the  form 


where  Q  is  the  stress  required  to  move  the  dislocation  from  its  original 
o 

lattice  site  to  the  point  of  neutral  equilibrium. 

Thus,  an  approach  to  the  problem  of  non-elastic  deformation  will  be 
developed  on  the  basis  that  the  mechanism  of  the  non- elastic  or  plastic  deforma¬ 
tion  is  the  formation  and  movement  of  dislocations. 

First,  a  detailed  mathematical  description  of  the  geometry  of  dislocations 
will  be  presented  and  the  relation  of  dislocations  to  internal  stress  and 
boiindary  conditions  will  be  demonstrated.  Then,  the  condition  for  movement 
of  the  dislocations  through  the  solid  will  be  formulated.  A  general  form 
of  the  equations  of  motion  based  on  this  mechanism  will  be  derived.  Repre¬ 
sentative  solutions  for  three  iniportant  problems  will  be  shown.  The  first 
solution  is  the  one- dimensional  rod.  The  second  solution  is  for  a  strip 
load  on  an  infinite  half- space.  The  physical  results  are  discussed  in  light 
of  the  features  of  the  deformation  sequence  and  observable  quantities. 


6 


The  third  prohlem  to  be  solved  is  the  defonnation  of  a  thin  cylindrical  shell. 
Approximate  equations  for  the  plastic  deformation  under  short  duration  loads 
will  be  set  up  and  solved. 
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II.  Ifathematlcal  Theory  of  Dislocations 


II. 1  Geometrical  Description  of  Dislocation 

The  deformation  of  an  imperfect  crystal  or  polycrystalline  medium  has 

been  recently  studied  utilizing  the  geometrical  descriptions  afforded  by 

(12)  (is) 

Riemannian  and  non-Riemannian  geometry.  The  work  of  Kondo'  ”  '  '  in  Japan, 

(14)  (2) 

Kroner'  in  Germany,  and  Bilhy  and  co-workers  in  the  United  Kingdom  have 
provided  a  great  Insight  into  the  theories  and  mechanisms  of  plastic  deforma¬ 
tion.  The  basic  mathematic  tools  that  these  investigators  have  utilized  is 
the  infinitesimal  parallel  displacement  transformation  which  forms  the  corner 

stone  of  the  general  theory  of  relativity  as  developed  by  Einstein  and 

(6)  (3) 

Kaufmann  and  Cartan  . 

In  this  approach  to  the  studies  of  a  plastically  deformed  crystal,  the 
stradned  crystal,  (  its  deformation  consists  of  elastic  strain  plus  added  dis¬ 
placement  due  to  dislocations )  or  the  xmstrained  crystal,  (vXrith  a  localized 
internal  strain  due  to  the  presence  of  dislocations.)  is  viewed,  with  respect 
to  a  reference  measuring  coordinate  system,  under  an  infinitesimal  parallel 
displacement.  This  transformation  may  be  shown  to  be  a  function  of  the  basic 
vectors  of  the  reference  system  times  a  connection  which  does  not  have 
tensor  properties.  If  the  metric  tensor  of  the  crystal,  whether  \mder  stress 
or  not,  is  synmetric  (indicating  only  elastic  deformation)  this  connection 
is  symmetric  and  shows  no  tensor  properties  whatsoever.  However,  if  the 
metric  tensor  which  describes  the  original  configuration  before  parallel  dis¬ 
placement  is  not  symmetric,  the  connection  is  not  symmetric.  Under  this 
condition  the  connection  may  be  uniquely  deconqMDsed  into  a  symmetric  part 
which  again  has  the  properties  of  a  connection  and  an  antisymmetric  i)art 
which  is  an  antisymmetric  tensor.  This  tensor,  which  is  called  the  torsion 
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becomes  an  in^Kirtant  quantity  in  describing  the 


(3) 

tensor  after  Cartan, ' 
nature  of  plastic  deformation. 

To  describe  the-  deformation  of  a  medium  with  dislocations,  the  connec¬ 
tion  and  torsion  tensor  must  be  evaluated. 

Consider  a  solid  body  which  is  defined  in  a  reference  Euclidian  three- 
space.  If  the  body  is  a  i)erfect  crystal,  l,e,,  no  inq)erfections  of  any 
natvire  are  present,  a  Cartesian  coordinate  system  may  be  attached  to  the 
medium.  Under  external  forces,  the  medium  goes  over  into  a  strained  configura¬ 
tion.  The  Cartesian  coordinate  system  maps  into  a  curvilinear  coordinate 
system.  If  the  basic  vectors  of  these  two  coordinate  systems  are  denoted  by 
for  the  Cartesian  system  and  e^  for  the  deformed  system,  metric  tensors 


for  these  coordinates  may  be  formed. 

hi  "  -i  *  (Il.l.l) 

‘  £ j  (II. 1.2) 

If  the  length  of  an  element  of  arc  in  the  original  and  deformed  systems 

ds^  "  (II.I'S) 

ds^  =  dx^  dx*^  (II.1.4) 

are  coiiq>ared; 

ds^  -  dsj  =  (g^j  dx^  dx^  (II.1.5) 


The  elastic  strain  tensor  is  defined  naturally  by 
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If,  however,  the  deformation  occurs  in  a  real  crystal  or  in  a  poly¬ 
crystalline  aggregate,  two  sets  of  basic  vectors  are  available.  The  first 

set  in  the  Cartesian  reference  frame  x..  The  second  set  of  basic  vectors 

—1 

is  given  by  the  basic  vectors  of  the  crystal  lattice,  denoted  by  e.^^  .  If 
there  are  no  imperfections  in  the  crystal  or  no  grain  boundaries,  these  two 
sets  of  basic  vectors  are  identical.  However,  if  the  clystal  is  not  perfect, 
either  initially  or  at  some  later  stage  of  the  deformation  process,  the 
unique  one-one  correspondence  is  no  longer  there  and  these  two  sets  are 
different.  It  becomes  apparent,  then,  that  the  description  of  the  deforma¬ 
tion  process  is  different  in  the  two  coordinate  systems.  If  the  deformation 
is  described  with  respect  to  the  Cartesian  basic  set,  the  total  macroscopic 
strain  will  be  shown.  If  the  second  set  is  used,  the  strain  that  is  shown 
would  be  just  the  elastic  distortion  of  the  crystal. 

Comparing  these  two,  the  difference  in  the  descriptions  should  give  a 
measure  of  deformation  introduced  by  non-elastic  means,  i.e.,  the  imperfections. 

This  difference  for  the  displacements  is  established  by  considering  an 
integral  expression  around  a  closed  path  in  the  real  crystal  and  the  corres¬ 
ponding  path  in  an  ideal  or  perfect  crystal. 

If  ^  represents  a  vector  from  x^  to  x^  +  dx^,  then  in  terms  of  the 
lattice  basis 


k 


e,  CO 

,1C 


(II. 1.6) 


where 


Se 


dx" 


dx" 


(II. 1.7) 
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is  the  Pfaffian.  The  expression  in  the  "brackets  may  or  may  not  be  integrable. 
The  change  in  the  lattice  basis  vector  is 


k 


iix.  Sk 


k  i 

where  the  F  is  the  connection  of  parallel  transformation  from  x 


.kv 


then 


^k 


Consider  the  expression 


If  the  expression 


a^M 


a% 


ax^Sx"^  Sx'^  ^x^  ^x  ^ 


0 


the  deformation  is  compatible,  i.e.,  there  are  no  imperfections  in  the 
crystal . 


to  e^^ 
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Now 


^  A*^  +  e 

■  axJ  J 


ik  a‘‘j  where  Sj 


-  ^  ^  -k 


hx 


Ik 


] 


If  the  vector  dM  is  integrable  around  the  closed  path 


a 


j 


If  the  vector  ^  Is  not  integrable  around  the  closed  path  in  the 
in5)erfect  crystal 


a^M 

axPix*^  ^ 


®k 


The  non-vanishing  connection  J’  [^X] 


(II. 1.12) 

-  (II.1.13) 


(IT.1.14) 


(II. 1.15) 


(II. 1.16) 


(II. 1.17) 


(II. 1.13) 
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The  component  of  normal  curvature  has  been  expressed  in  the  original  Euclidian 
coordinates,  A  similar  expression  may  be  derived  for  the  normal  curvature  of 
the  nature  or  lattice  basis. 


h  = 


V 


(>  Ji  p^  pP  +  A  aP 

^  V  IJL\  ^  ^  JUp^  pX 


(n.1.19) 


Coii5)Uting 

SxtJax^i 


ax^axM 


(II. 1.20) 


Thus  the  antisymmetric  part  of  the  connection  is  not  zero.  This  portion 
has  tensor  properties  iinlike  its  symmetric  counterjjart .  Denote  the  tensor 
(the  torsion  tensor  of  Cartan), 

pk  =  gk 

[,pX]  ^  ^  pX 

Since  the  basic  mechanism  of  plastic  deformation  in  the  theory  of  Taylor 
involves  the  effects  of  a  dislocation  and/or  the  formation  and  transport  of 
these  dislocations  through  the  medi\im,  the  question  may  be  asked,  what  does 
this  three-index  tensor  have  to  do  with  the  deformation?  Obviously,  the 
formation  and  movement  of  the  dislocation  is  not  included  in  the  formulation. 
Thus,  this  description  must  describe  the  dislocation  itself. 

The  topology  of  the  crystal  is  indicated  by  the  difference  in  the  outward 
normal  vector  of  the  perfect  crystal  circuit  and  the  outward  normal  vector  of 
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I 


J 

i 

the  corresponding  vector  in  the  imperfect  crystal.  This  is  the  Burger's  vector, 
which  is  defined  as  the  difference  between  the  outward  normal  vector  of  two 
circuits  in  the  material  (See  Figure  2),  the  first  circuit  being  a  closed  path 
in  the  iiiQ)erfect  material  and  the  second  the  one-one  corresponding  path  in  the 
ideal  crystal. 

Consider  these  two  circuits  in  the  reference  plaixe  and  the  mapped  plane.  (Fig. 
The  transformation  from  imperfect  to  perfect  crystals  consists  of  the 
translation 

=  2  032^  (II. 1.21) 

plus  the  rotation 

=  "Vx  (II.1.22) 

Then  a  vector  ^  at  x  of  the  original  frames  goes  into  T  after  one  circuit 
in  the  perfect  frame 

I  *•  i  =  R  VMX  i  CDi  0)2  (II. 1.23) 


The  Burger's  vector  is 


b  = 


(II. 1.24) 


-■k 


.k 


The  quantities  ^ 


,  identify  the  torsion  tensor  and 
curvature,  which  define  the  imperfection. 


If  these  quantities  are  examined  separately  the  physical  interpretation 

k 

is  obvious.  Assuming  that  the  metric  is  everywhere  compatible,  i.e.,  R  =0- 


l4 


Burger's  Circuits  Around  a  Dislocation 


Figure  2 
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V  a 
a>i  (X)z 


(II.I.25) 


k 

Only  two  essentially  different  components,  S  ,  of  the  torsion  tensor  are 

fj,  {y 

/I  2  ^ 

noticed  if  the  mediiim  may  be  regarded  as  isotropic.  (S  s  ®  12^  (See  Fig,  4) 

k  ./ 

The  torsion  describes  the  local  behavior.  On  the  other  hand,  assume  R  0 

Then  four  cases  arise  (see  Figure  5)* 


1. 

^121" 

122 

2. 

r\21^' 

«  122 

3, 

pi 

^123 

k. 

^  123 

From  the  expression  for  the  Burger's  vector  the  first  corresponds  to  the 
slip  of  an  edge  dislocation  not  parallel  to  the  edge,  the  third  corresponds  to 
slip  of  an  edge  dislocation  parallel  to  the  edge  and  the  second  and  fourth 
correspond  to  the  slip  of  a  screw  dislocation  not  parallel  to  the  edge  of  the 
screw,  respectively. 

The  curvature  term  denotes  in5)erfection8  in  the  large.  (Figure  5)» 


11,2  The  Dislocation  Density 

From  the  foregoing  discussion,  the  curvature  of  the  system  denotes  imper¬ 
fections  in  the  large.  The  torsion  tensor  denotes  inQjerfections  in  the  small 
and  actually  ref ex's  to  the  dislocation,  which  is  the  basis  of  the  theory  of 
plasticity.  Noting  that  descriptions  of  plastic  deformation  usually  do  not 
involve  descriptions  in  the  large  and  that  the  material  that  is  being  deformed 
plastically  maintains  macroscopic  continuity,  the  effects  of  a  non- zero 
curvature  will  be  ignored.  Attention  will  be  focused  on  the  torsion  tensor. 
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If  the  discrete  lattice  basis  of  the  real  crystal  is  replaced  by  a  con¬ 
tinuous  lattice  which  has  the  same  metric  properties  as  the  lattice  basis, 
certain  additional  relations  may  be  obtained  by  considering  the  differential 
relations  concerning  the  Burger's  vector. 

The  Burger's  vector  is 

^  (II.2.1) 

The  differential  expression  for  the  Burger's  vector  is 

^  (II. 2. 2) 

(2)  .  (ik) 

Bilby  and  Kroner  have  demonstrated  that  this  tensor  may  be  written 
as  the  curl  of  second  order  tensor  if  the  torsion  tensor  vanishes  at  infinity. 
Introduce  the  term  defined  by 

^  (11,2.2) 

where  is  the  Eddington  nxunber 


ppy 


even' 

odd 


permutation  of  1,2,3 


otherwise 


Then 

e,  db^  =  L^^e  de 

— k  -p  — k 


(II. 2. 4) 
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k  p 

The  quantity  L  ^  is  a  second  order  tensor  relating  the  Burger's  vector 
to  the  natural  "basis  and  is  called  the  dislocation  density.  An  analogous 
definition  may  "be  made  with  reference  to  the  Cartesian  reference  "basis.  This 
quantity  corresponds  to  the  specification  of  an  ar"bitrary  distri"bution  of  dis¬ 
locations  in  the  medium. 

The  simple  physical  interpretation  of  the  dislocation  density  is  that  it 
represents  the  Burger's  vector  in  the  k  direction  from  a  single  dislocation 
which  lies  in  the  p  direction. 

The  dislocation  density  is  treated  as  a  continuous  function,  defined  by 
the  limiting  process  above.  Another  interpretation  which  will  have  use  in 
the  application  may  be  obtained  if  the  dislocation  density  tensor  is  obtained 
by  forming  the  torsion  tensor  around  each  discrete  dislocation  and  then 
foimiing  the  dislocation  density  by  the  same  limit  process.  This  results  in 
a  sequence  of  non-zero  values  at  isolated  points.  Then  the  complete  dis¬ 
location  density  tensor  may  be  formed  by  taking  the  average  value. 


i  1 

.kp  ^  ^  ^(v)  ^  (v)) 

Z  cos  (x^  - 


(11.2.5) 


k  jO 

where  the  dislocation  density  associated  with  a  discrete  disloca¬ 

tion  Vf  k^  is  the  reciprocal  of  the  distance  between  dislocations  in  the  1 

direction,  is  the  reference  coordinate,  and  is  the  coordinate  of 

the  dislocation. 

Passing  to  the  limit  of  a  continuously  dislocated  crystal 
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(II. 2. 6) 


kp  -I 

'  o> 

/  cos  (k^x^)  dk^ 


The  operational  significance  of  a  derivative  of  some  quantity  containing 
the  dislocation  density  in  the  reference  plane  now  may  he  interpreted  as 
division  hy  the  average  distance  between  dislocations. 


II. 3  Stress  and  Dislocations 

The  next  step  in  the  development  of  dislocation  theory  of  plastic  deforma¬ 
tion  is  to  derive  the  relation  between  stress,  elastic  strain,  the  dislocations, 
and  their  relationship  to  boundary  tractions. 

The  basic  approach  to  the  stress-dislocation  relationship  is  made  through 
the  following  procedure.  Consider  a  perfect  crj-^stal  with  no  surface  tractions. 
There  is  no  stress  or  strain  in  the  interior  of  this  crystal  by  definition. 

Now  insert  a  dislocation  into  this  crystal,  keepiiag  the  surface  tractions 
zero.  Now  there  is  a  definite  non-zero  strain  in  the  sjjecimen.  An  inter¬ 
nal  stress  which  corresponds  to  the  strain  is  generated.  This  is  the 
internal  stress  due  to  the  dislocation. 

To  demonstrate  this  mathematically,  the  deformation  must  be  viewed 
from  both  the  reference  or  macroscopic  coordinate  system  and  from  the  lattice 
of  natural  coordinate  system.  Thus  the  concept  of  shape  or  total  deforma¬ 
tion  and  lattice  defcymiation  is  introduced. 

From  the  definition  of  strain,  the  strain  may  be  written  in  terms  of 
the  connection  and  the  basis  vectors  of  the  reference  coordinate  system 


(II. 3.1) 
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From  the  above  demonstration,  the  connection  consists  of  a  symmetric  part 
and  an  antisymmetric  part,  thus  giving  rise  to  pure  strain  and  rotation. 

The  strain  in  the  lattice  coordinate  system  has  a  similar  development 


(n.3.2) 


If  this  lattice  strain  is  written  in  the  reference  coordinate  system 


^  (L)  _  Ml  k 
Ij  '  ^  Jk  ^ 


(II. 3. 3) 


where  the  connection  is  given  by  transformation  of  coordinates. 

Since  the  total  deformation  consists  of  the  lattice  strain  and  disloca¬ 
tion,  the  deformation  may  be  written  as  a  sum 


+  T 


(D) 


(II. 3. 


Introducing  a  similar  connection  for  the  deformation  from  the  dis¬ 
locations,  it  is  readily  seen  that  the  connections  are  additive.  The 
symmetric  part  of  the  connection  gives  rise  to  a  strain  and  the  anti¬ 
symmetric  part  to  a  rigid  rotation.  Since  the  con5)onents  of  this  decom¬ 
position  are  independent  of  the  other  component,  it  is  evident  that  the 
total  strain  is  a  sum  of  the  lattice  strain  and  the  strain  is  Introduced  by 
the  dislocation. 

To  incorporate  a  discussion  of  stress  into  dislocation  theory,  it  is 
reasonable  to  assume  that  the  stress  is  a  function  of  the  lattice  strain 


a  -  M 

®  ki 


(II. 3. 5) 
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ki 

where  M  are  ela-stic  coefficients, 

1 J 

The  lattice  deformation  is  given  hy 


T  — 


For  differential  Increments  of  deformation 


(11.3.6) 


dx'' 

If  the  strains  are  small,  the  curvature  is  given  by 


^  Jki  ^k^  Ji  "  jk 


(11-3.7) 


(11.3,8) 


Associated  with  any  Burger's  circuit  there  is  a  change  in  lattice  deformation 

k£ 

where  X  is  the  differential  of  area. 

If  the  lattice  strain  may  be  written  as  a  fvmction  of  position,  the 
ctirvature  tensor  may  be  contracted 


(II. 3. 10) 


where 


r  =  c 

a]3v  ®aX 


(11,3.11) 


and  e  is  the  Eddington  symbol.  Since  the  lattice  strain  and  rotations  are 
fimctions  of  the  connections  and  since  the  connection  may  be  written  as  a 
function  of  a  second  order  tensor 

k 


(II. 3.12) 
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(II. 3.13) 


J  ]k 


dx' 


k 


1 


?7  64  L  —  Lj  ,.  +  e 


ij  '  ""la^  ^ 


(11.3.14) 


The  contracted  curvature  tensor  is 


^Ij  ”  ^ja|3  ^176  ^/3  "^ja|8 


(II. 3. 15) 


If  both  the  lattice  strain  and  lattice  rotations  are  functions  of  posi¬ 
tion,  the  lattice  deformation  is  compatible  in  the  large  and  the  curvature  is 
zero.  If  the  linear  elastic  stress- strain  law  is  valid,  a  compatibility  equa¬ 
tion  between  stress  and  dislocation  density  follows 


'la/3  J75  ^^7 


-e 


(11.3.16) 


Solutions  to  this  problem  for  arbitrary  surface  tractions  and  boundary- 
conditions  are  given  by  the  following  procedure.  For  this  it  is  assumed  that 
the  right  hand  side  is  known  as  a  function  of  position.  Denoting  the  right 
hand  side  by  the  solution  for  the  strain  is  given  by^^^ 

(r  -  r') 


IJ  Wn  V 


dV 


(II. 3. 17) 


Then  the  elastic  stresses  and  strains  for  the  perfectly  elastic  body  of  the 
same  shape  subject  to  the  same  boundary  conditions  are  determined.  Denoting 
this  solution  by  "the  f\all  solution  is 


®ij  ®'ij  ®ij 


(L) 


(11.3.18) 
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II. 4  Stresses  and  Moviiig  Dislocations 


The  foregoing  discussion  lays  a  detailed  foundation  for  the  study  of 
plastic  deformation.  The  introduction  of  the  torsion  tensor  as  a  measure  of 
the  internal  forces  generated  hy  a  dislocation  is  a  key  point.  This  shows  that 
the  existence  of  dislocations  in  otherwise  perfect  crystals  is  possible  and 
indeed  quantitatively  compatible  with  any  discussion  of  inelastic  deformation. 

The  foundation  has  been  restricted  to  static  deformation.  As  mentioned 
in  the  introduction,  the  basis  of  Taylor's  original  theory  was  that  plastic 
deformation  is  basically  caused  by  the  movement  of  dislocations  through  the 
medium  imder  stress. 

A  natural  generalization  to  this  development  then  is  to  introduce  a 
quasi-torsion  tensor  which  Incorporates  both  the  internal  stress  around  the 


dislocation  and  the  movement  of  the  dislocation. 


Consider  the  following  three-index  quantity. 


■  MV 


(11.4,1) 


Let  us  examine  the  situation  if  several  dislocations  pass  through  a 
unit  voliime.  Initially  the  volume  contains  one  dislocation,  Define  the 
Burger's  circuit  as  before.  The  tensor  is  Just  the  torsion  tensor 

of  the  previous  section. 


Jc 


(II. 4. 2) 


Now  allow  the  dislocation  to  move  through  and  out  of  the  voliime  element, 
and  let  another  dislocation  enter  into  this  volume  and  assume  the  same 
position  (see  Figure  6). 
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lAirgftr'r;  CirculLs  for  Moving  Dlsl 
Figure  6 
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The  original  Burger's  circuit  has  now  become  two  separate  non- closed 
paths.  The  Burger's  vector  of  the  dislocation  has  remained  unchanged. 
However  the  displacement  of  the  entire  upper  half  of  the  volume  is  one 
lattice  unit.  Denote  this  displacement  by  a  vector  f  . 

The  tensor  H  may  be  written  then  in  terms  of  this  vector,  which 
has  unit  length  per  dislocation  moved  through  the  volume  and  the  number  of 
dislocations. 


ok  ^  Sr 


(II. “1.3) 


where  n  is  the  niimber  of  dislocations  that  have  moved  in  the  vth  direction. 
V 


This  quantity,  r) 


hf 


k 


,  is  the  flux  of  dislocation  in  the  v-th  direction 


times  the  strain  caused  by  the  elementary  movement  of  a  dislocation  in  the 
p-*  th  direction. 

It  is  readily  seen  that,  if  many  dislocations  move  through  the  body,  the 
displacement  arising  from  the  flux  term  will  be  greater  than  the  local  strain 
induced  by  the  presence  of  the  dislocation. 

Thus  an  approximation  for  plastic  flow  presents  itself.  The  local 
strain  field  due  to  the  presence  of  the  dislocation  is  ignored  and  only  the 
terms  arising  from  the  flux  of  the  dislocation  are  retained. 

If  the  general  flux  term  is  treated  as  a  generalized  flux  in  the  Onsager 
sense,  equations  may  be  established  for  combined  elastic-plastic  deformation. 
The  force  that  corresponds  to  this  dislocation  flux  is  taken  'by  analogy  to 
diffusion  of  point  vacancies  as  gradient  of  the  dislocation  density. 

The  crux  of  the  actual  dynamic  plastic  mechanism  is  that  the  disloca¬ 
tions  move  under  stress,  If  the  dislocation  is  viewed  end  on,  the  movement 
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of  the  dislocation  is  analogous  to  diffusion  under  stress  of  a  vacancy.  It 
should  then  be  seen  that  the  movement  of  a  dislocation  shoxxld  be  dependent 
on  the  stress,  the  temperature,  and  the  concentration  gradients  of  disloca¬ 
tion.  Somewhat  heuristically,  the  movement  of  the  dislocation  will  be  treated 
similarly  as  a  diffusing  vacancy.  The  dislocation  will  move  if  the  vibration 
of  the  atom  on  the  edge  of  the  dislocation  is  large  enough  in  the  direction 
of  the  lattice  imperfection  caused  by  the  dislocation.  In  addition,  all  of 
the  atoms  along  the  dislocation  must  satisfy  this  requirement  simultaneously. 
The  probability  is  affected  by  the  neighboring  atoms  which  must  move  away 
sufficiently  to  allow  passage  of  the  atom.  The  probability  that  a  dislocation 
will  move  under  an  applied  stress  field,  then,  is  a  compound  probability,  com¬ 
posed  of  the  two  factors; 

1,  The  probability  that  one  atom  on  the  dislocation  will  move. 

2.  The  probability  that  two  adjacent  atoms  will  have  the  same 
probability  for  a  jump  at  the  same  time. 

D\ 

The  diffusion  coefficient  of  the  dislocation  may  be  written' 

D=ra^  (II.J+.4) 

where r  is  the  probability  of  a  jump,  and  a  is  the  jump  distance 


-  -U/kT-tr^  Tj/kT 
r  =  V  e  &  IT  e  V 


^oi  ^ol 


(II. '+.5) 


where T  is  a  weighted  thermal  frequency,  is  the  activation  energy  for  the 
atom  which  is  going  to  jiimp,  U.  are  the  activation  energy  contributions  caused 
by  the  neighboring  atoms,  is  the  pair  interaction  potential,  and  is 

the  probability  that  the  1-th  atom  will  be  in  the  same  jump  configuration. 

The  barrier  is  greatly  reduced  by  a  shearing  stress.  The  probability  of 
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the  neighboring  particle  being  in  jump  position  increases  as  the  shearing  stress 


lowers  its  activation  barrier.  Assuming  that  the  are  not  affected  by  the 
shear,  an  approximate  expression  for  the  diffusion  coefficient  may  be  written  if 
the  barrier  energy  is  a  linear  function  of  distance  and  is  expanded  in  a 
Taylor  series  in  the  shearing  stress.  Retaining  only  the  first  term  for 


-U 

-2+2^ 
T  kT  T 

"o'  ° 


(II.1I.6 


where 


D  =  a^jr e  i 


-U 


g. 


(iiA.t) 


o  ■  ^  '^oi 

T  is  the  shear  stress,  is  the  shear  stress  corresponding  to  a  strain  which 
moves  the  atom  to  the  top  of  the  barrier. 


Since 


- 

-TT  „  I?!’ 


r  ^  vl  e  -  g  —  e 
i  .  o 


T  „  T#  ^  TT 


(II.4.8) 


=  ^  Sol  77  ® 


-(Uo+SOl)  +  ^ 

-  o 


which  is  the  probability,  can  never  be  greater  than  one,  the  diffusion  coefficient 
has  an  upper  value  D  =  when  the  shear  stress  is  greater  than  . 

The  effect  of  hydrostatic  stress  or  pressure  on  the  yield  strength  of  the 
material  enters  thii*  ^expression  by  changing  the  magnitude  of  the  potentisd 
barrier.  To  a  first  approximation. 

U  (01+  02+  03)  =  +  ^3) 

^(cyi.  +0*2  ^^3/ 
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where  CTi,  (Ta  i  are  the  principal  stresses  and  aj.  +CT2  +CT3  is  three  times 
the  hydrostatic  pressure. 

The  exact  form  of  this  change  may  he  quite  involved  and  is  not  considered 

here. 

A  great  ntimber  of  dislocations  are  contained  in  a  typical  element  of  voltime. 
The  orientation  of  these  dislocations  for  motion  may  he  considered  as  randomly 
oriented.  If  this  is  the  case,  the  macroscopic  or  effective  diffusion  coefficient 
will  he  different  than  the  diffusion  coefficient  for  a  siiigle  dislocation. 

Assume  an  element  of  volvune  with  applied  stresses.  Orient  this  volume  so  that 
the  maximum  shearing  stress  is  along  the  X-axis.  Denote  the  orientation  of 
the  dislocation  movement  line  hy  6  measured  from  the  X-axis.  (Figure  7) 


^wax 


Since  the  shearing  stress  is  zero  for  0  =  45'*,  there  will  he  a  restricted 
zone  of  orientation  of  dislocations  that  will  move  rnider  this  stress,  say  up  to 
Tl/,Tf/  <  0,  Then  the  macroscopic  diffusion  coefficient  is  given  hy  averaging  the 
component  of  D  along  the  X-axis  over  all  possible  orientations.  Then 

D  =  D  (11.4.9) 

o  f 
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where  D  is  the  macroscopic  diffusion  coefficient,  and  is  the  single  dislocation 

diffusion.  For  45”  the  value  of  is  between  0.9  and  1.0.  Thus,  for 

Tj/ 

the  purposes  of  this  paper  D  will  be  taken  equal  to  D^, 

The  equation  of  motion  for  the  elastic  plastic  medium  may  be  written 


Sv.  So. . 

p _ 1  ^ 


(II. 4. 10) 


where  is  the  particle  velocity  in  the  i-th  direction. 

The  equations  of  motion  for  an  arbitrary  body  referred  to  a  reference  Cartesian 
system  may  be  written 


(II. 4. 11) 


where  a  are  the  contravariant  con^jonents  of  the  stress  tensor  given  in  the  local 
metric;  ij  refers  to  total  covariant  differentiation,  p  is  the  density  and  is 


the  total  or  convective  time  derivative.  The  total  covariant  differentiation  is 
defined  as 


nin 


(11.4.12) 


where  the  comma  denotes  partial  derivatives  andP  is  the  convexion  of  the 

mn 

metric  space.  In  the  case  of  dislocations,  the  metric  tensor  that  describes  the 
space  is  given  by  the  lattice  unit  vectors.  In  the  case  where  no  dislocations 
are  present,  the  lattice  vectors  coincide  with  the  Cartesian  reference  system 
and  the  connexion  is  identically  zero.  If  one  or  more  dislocations  are  present, 
the  lattice  is  distorted  and  the  connexion  is  given  by  the  E\iler-Schouten  tensor, 


i  i 

r  =  S  If  a  dislocation  has  passed  through  an  elementary  volume  in  the 

original  crystal  lattice,  the  original  coordinates  are  displaced  by  an  integral 
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nimiber  of  con5)lete  lattice  units.  This  displacement  is  described  by 

pi  =  N  ^ 
mn  n  ^  m 


(11.4.15) 


where  i^  the  number  of  dislocations  that  have  passed  through  the  elementary 

voltune  in  the  m-th  direction  and -  is  the  projection  of  this  displacement  in 

the  i-th  direction.  If  no  dislocations  are  present,  the  lattice  metric  coincides 
with  Cartesian  reference  metric  and  the  connexion  is  zero.  If  a  dislocation  is 
present,  the  connexion  is  given  by  the  Euler-Schouten  tensor.  If  a  dislocation  has 
passed  throi:igh  the  region  \mder  some  force  and  is  no  longer  present,  a  residual 
three  indiced  tensor  demonstrates  that  a  slip  has  occurred.  The  resulting 

connexion  is  the  siim  of  these  two  connexions 

(II. 4.14) 

The  amovmt  of  deformation  due  to  the  presence  of  a  dislocation  is  small  com¬ 


Tim  mn  mn 


pared  to  the  translatory  movement  caused  by  the  passage  of  dislocations. 

Now  examine  the  contravariant  derivative  of  the  stress  tensor 

_  J-J  4.  Tpi  (11.4.15) 


a-  =  a''  + 

~  ,j  mn 


If  there  are  no  applied  surface  tractions  the  stress  at  any  interior  point  of 
the  body  is  due  to  the  elastic  strain  surrounding  the  dislocations.  If  this  stress 
state  is  vinstable  thermodynamically,  the  dislocations  will  move.  If  the  surface 
tractions  are  zero 


and 


=  0 


=  -f" 

,j  mn 


(II.4.16) 


(11.4.17) 
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Thus  a  stress  gradient  on  the  Interior  of  the  body  exists  due  to  the  moving 
dislocations.  Since  Hooke's  law  for  elastic  strains  has  been  assumed,  addition 
of  STirface  tractions  to  the  body  adds  stresses  linearly  to  those  in  the  body  due 

to  dislocations.  Define  now  a  stress  tensor 

„mn  mn 
b  =  a 

which  denotes  the  stress  due  to  the  dislocations  only,  then  the  equation  of 
motion  may  be  written 


,  av^ 

=  p 


(II. 4.18) 


ax^  ax""  ■  at 

The  other  eqiiations  necessary  to  complete  the  description  of  the  deformation 
process  are  derived  from  energy  and  entropy  considerations.  Physically,  it  may 
be  expected  that  the  equation  for  motion  of  dislocations  through  a  solid  body  will 
obey  a  diffusion  equation.  To  investigate  the  transport  properties  of  the  disloc¬ 
ations,  we  note  that  the  change  in  energy  during  this  transport  is  proportional  to 
the  time'  rate  of  change  of  the  total  strain.  Similarly  the  energy  change  may  be 
written  as  a  change  in  the  connexion  with  position.  Equating  these  two  quantities 
should  lead  to  the  subsidiary  equations  necessary  to  coii5)lete  the  description  of 
the  motion. 

The  energy  equation  may  be  derived  as  follows:  We  consider  the  following 
differential  expansion  for  the  internal  energy  in  terms  of  the  lattice  strain  of 
the  dislocated  crystal. 


-  "ij  *«  * 


(11.4.19) 


where  E  is  the  internal  energy  per  \init  mass,  d  is  the  stress,  ^  is  the 
lattice  strain  and  S  is  the  entropy  per  unit  mass.  The  strain  may  be . deconposed 
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into  parts  arising  from  surface  tractions  and  the  dislocations. 

The  lattice  strains  due  to  the  movement  of  dislocation  may  be  written  as 


de,  ,  =  dx^ 

ij  jk - 


where  ^  is  the  basis  vector. 

If  the  strains  are  small,  the  strain  may  be  written  as 

de.  .  =  ^  F^  .  -  ^  F^ .  .  dx^  dx^  =  ,  dx^  dx^ 

ij  1  jk  k  ij  jkl 


(II. 4. 20) 


Since  the  medivun  luider  consideration  is  isotropic,  we  may  contract  the 
curvature  term. 

k  ,  k 


de.  .  =  R.  ,  dx  dx  = -e,  e.  ^ 

ij  1  jkk  j  mn  '^isp  °m  On  : 


isp 


np 


^  A  L-L  ''  (11.4.22) 

“  J  mn  ^m  '  ^n  in'” 

where  is  the  Eddington  permutation  symbol  and  is  the  dislocation  density. 

If  the  reversible  part  of  the  strain  energy  is  subtracted  from  the  energy  equation 

(11.4.23) 

Since  the  process  is  diffusive  in  nature,  we  postiilate 


TdS.  =  e  .  a  (a.  L  -  L  ) 
irr  jmn  m  in  in 


a^^inr* 

DT  -  =  T 

at 


a  irr 


where  B  is  the  stress  dependent  diffusion  coefficient  for  dislocation  movement. 
Inserting  this  expression  for  entropy  and  contracting  indices 


The  dislocation  density  is  proportional  to  the  plastic  strain  as  a  first  approx¬ 
imation  and  thus 


D 


e.  . 


(II.'+.25) 


Resolving  the  strain  into  hydrostatic  and  deviator  con^jonents 


e. .  =  6 .  .  0 

ij  ij 


+  d 


ij 


and  adding  and  subtracting 


Di 


3x^3x^ 
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(IIA.26) 


(11.4.27) 


^  ”  dt 


(II. 4. 28) 


The  equations  of  Motion  are 


Sa. .  ^S._  3v. 

+  p_l 

5x  4  dx_ 


m 


(11.4.29) 


72  ^  ^ 


(11.4.30) 


Da  d. 


"  at 


(Il.i^.31) 
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III.  Ing)act  Plastic  Deformation  of  a  One  Dimensional  Rod 


The  first  problem  to  be  solved  for  dynamic  plastic  deformation  is  the  case 
of  a  one  dimensional  rod  in5)acted  at  the  end.  The  rod  is  considered  as  a  con¬ 
tinuum  extending  from  0  to  infinity  in  the  x  direction.  A  load  is  applied  at 
the  end  x  =  0  and  is  given  by  (t). 

The  equations  for  one  dimensional  dynamic  plastic  deformation,  where  the 
elastic  strain  is  considered  to  be  linear,  are 


jij  ^ 

'  ax  ”  ^  at^ 


(m.i) 


afu  ^  a^u 
ax®  axat 


(III.2) 


where  M  is  an  elastic  constant 
p  is  the  density 
u  is  the  displacement 
D  is  a  diffusion  coefficient  and 

a  is  the  plastic  driving  force  or  the  stress  due  to  the  dislocations 
defined  by  expanding  this  as  a  Taylor  series  in  the  total  strain 


0 


ax^  ax 


(111.3) 


and  a  and  jS  are  the  constants  relating  the  stress  and  dislocation  density. 
a,/3  and  D  are  dependent  upon  material,  temperature  and  stress  configuration. 


D  =  Dq  7"  ( 

o 


T  ^  T  = 
O 


t.AH 

I  • 
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a=  a 

o 

D  =  D 

o 


T  >T 
O 


T  ^ 
c  RT 


(III. 5) 


where  T is  the  maximum  elastic  shear  stress.  In  this  problem  which  is  a  case 
of  plane  strain,  the  lateral  stresses  are  prescribed  and  thus  the  maximum  shear 
stress  in  these  expressions  may  be  replaced  by  an  equivalent  expression  in  the 
elastic  strain. 


The  first  of  the  equations  represents  the  balance  of  the  forces  due  to 
elastic  strain  particle  acceleration  and  the  dislocations,  and  the  second  equation 
represents  the  dissipating  of  energy  by  motion  of  the  dislocations  through  the 
crystal  lattice.  The  energy  that  is  dissipated  shows  up  as  thermal  energy  in  the 
crystal.  The  internal  stress  or  body  force  due  to  the  presence  of  the  dislocations 
is  called  the  plastic  driving  force. 

The  diffusion  coefficient  may  more  accurately  be  called  a  kinematic  viscosity 
but  this  distinction  is  semantic. 

Thus  the  first  equation  may  be  written 


M  g  ^  ^^u  _  ^^u 


(III. 6) 


a 


Multiply  the  second  equation  by  X  =  — 
Xd  ^  =  x& 


ax^ 


at^^ 


(III.T) 


and  subtract  from  the  first  equation. 
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I 


Then 


u 


^  -X 


p  ax' 


at' 


axat 


(III. 8) 


or 


a^u  _ 

C2 

.  fx 

a^u 

at^ 

ax^j 

1  axat 

pD 

ax^ 

i 

(III.9) 


The  first  term  in  this  equation  corresponds  to  the  elastic  wave  equation. 
By  itself  an  undanq)ed  linear  elastic  wave  would  propagate  down  the  rod. 

The  second  term  is  the  contribution  from  the  dislocation  movement. 

Rewrite  the  second  term  as 


ax 


au  ^  ^ ^ 

at  pDx  ax 


(III. 10) 


Since  X=  ■2^  the  term  becomes 


X ^  ^ 

ax  at  a  ax 


(III. 11) 


Now  the  differential  operator  in  front  of  the  brackets  may  be  interpreted 
as  division  by  the  average  distance  between  dislocations, 
a  A  a_ 
d 


ax 


(III. 12) 


o  o 


A 

d 


^  +  a  D  ^ 


at 


ax 


(III. 13) 
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However 


m 

a 


may  be  interpreted  as  the  velocity  of  a  moving  dislocation. 


A  (  ^  +  V  ^ 
d  ^  St  Sx 


(III. Ilf) 


The  differential  equation  of  motion  may  be  written 


Su  ^  ^ 

St  “  Sx 


=  0 


(111.15) 


Consider  the  problem  imder  the  restriction  c>v>-CjV>  o.  This  is  a 
problem  which  has  been  investigated  fully  by  Whitham,  The  boundary  condition 
is  given 

P  =  J^Q(t)  t  =  0  X  =  0 

.^=u  =  0  t=0  x>0 

St 

The  second  boixndary  condition,  which  is  normally  the  displacement  at  in¬ 
finity  for  the  elastic  wave  propagation  in  the  rod,  is  modified  here.  Since 
the  plastic  deformation  is  caused  by  freeing  of  dislocations  under  stress  or 
strain,  the  other  boundary  condition  must  be  given  as  a  function  of  the  incident 
stress  wave,  This  condition  will  be  explicitly  evaluated  later.  It  sufficies 
here  to  state  that  the  elastic  wave  from  infinity  is  identically  zero.  The 

plastic  wave  is  caused  physically  by  the  movement  of  dislocations  freed  by  the 

state  of  elastic  stress.  This  second  boundary  condition  thus  is  an  internal 
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boundary  condition  and  is  treated  as  a  free  boundary  on  the  elastic  wave  front. 

The  first  boundary  condition  is  expressed’ as  a  pressure  or  stress.  This 
affects  only  the  elastic  strain,  not  the  total  strain.  The  total  strain  is  deter¬ 
mined  as  the  sum  of  elastic  displacement  and  the  flow  due  to  dislocation  movement. 

If  the  pressure  is  held  constant  the  elastic  strain  is  constant.  When  a  dislocation, 
however,  comes  to  the  surface,  an  elastic  rebound  may  be  expected,  and  this  re¬ 
bound  adds  an  elastic  strain  in  addition  to  the  strain  already  present  due  to 
the  pressure.  Because  of  the  continuum  notation  for  dislocations  this  rebound 
is  very  small  for  each  individual  dislocation,  and  is  continuously  spread  out 
on  the  bo\indary.  It  is  therefore  indistinugishable  from  the  already  present 
elastic  strain  and  can  be  safely  ignored.  This  problem  will  crop  up  again  in 
the  two  dimensional  problem.  There  it  will  be  stated  that  the  stress  due  to 
dislocations  is  zero  on  the  surface. 

The  displacement  at  x  =  0  is  composed  of  two  parts,  the  first  is  the  elastic 
strain  and  the  other  is  the  displacement  caused  by  the  movement  of  dislocation. 

Thus  the  boundary  condition  may  be  written  as  a  displacement 


u  =  u  (t)  =  u^^^^  (t)  (t)  t>0  x  =  0 

Denote  the  one  sided  Laplace  transform  of  the  displacement  by 


u(x,t)  dt 

and  the  inverse  Laplace  transform  by 


u(x, s)=s/  e 
o 


(III.I6) 


1  A. 

U  (x,t)  ="2^  ^Br  ^ (III. 17) 

where  Br  denotes  the  Bromwich  contour  in  the  complex  s-plane,  the  s  in  the 
denominator  is  the  resiilt  of  the  Heaviside  step  function  in  the  original  two 


sided  Laplace  Integral. 


4o 


The  differential  equation  in  the  transform  variable  may  be  written 


3“^  -  c‘ 


il 


5x 


u  =  0 


(III  i8) 


The  general  solution  of  the  transformed  displacement  is 


-  -fn(s)x 

u  (xss)  =  (s)  e 


+  (s)  e 


(EIJ.  ].9) 


where  f^  (s),  f^  (s)  are  solutions  of  the  auxiliary  equation . 


c  (f(s))  +  -  “3 


f(s)  -s(6  +  |)  =  0 


(HI.  20) 


For  large  s 


(III.  21) 


so  that  for  large  distances  from  the  end  of  the  rod  the  first  term  represents 
and  elastic  wave  propagating  down  the  rod.  The  second  term  represents  an  in¬ 
coming  elastic  wave  from  Infinity. 

For  small  s 


^ 


(111.22) 


The  first  of  these  corresponds  to  a  boundary  layer  at  x  =  0.  This  boundary 
layer  expresses  the  difference  in  the  actual  displacement  and  the  elastic 
component.  This  boundary  layer  term  is  necessary  to  fully  complete  the  boundary 
conditions  stated. 

The  second  term  corresponds  to  a  "wave"  being  propagated  backwards  with  the 
velocity  of  a  moving  dislocation. 


The  hovmdary  conditions  in  terms  of  displacements  may  he  further  clarified. 
The  presence  of  a  boundary  layer  at  x  =  0  implies  that  two  boundary  conditions 
are  necessary  to  specify  the  solution  to  the  problem.  The  first  of  these  boundary 
conditions  is  given  by  the  elastic  displacement  which  is  a  function  of  the  stress 
applied  at  the  surface.  The  difference  between  the  actual  displacement  and  the 
elastic  displacement  is  given  by  the  boundary  layer. 

This  added  displacement  arises  from  the  movement  of  dislocations  under  stress 
to  the  ends  of  the  rod;  .  These  dislocations  must  have  a  starting  point  for  motion. 
This  starting  point  can  be  none  other  than  the  front  of  the  elastic  wave,  for 
the  rod  does  not  know  that  the  dislocations  should  move  until  the  stress  wave 
arrives. 

Thus  the  second  boundary  condition  at  the  surface  may  be  replaced  by  a 
boundary  condition  along  x  -  ct, 

This  boTondary  condition  is  just  the  displacement  at  the  point  (x, ct).  This 
means  that  the  state  of  elastic  strain  on  the  interior  of  the  rod  determines  the 
movement  of  a  dislocation  at  that  point, 

It  is  seen  then  that  the  boimdary  condition  that  is  applied  at  the  surface 
is  just  the  elastic  displacement.  The  plastic  conqponent  of  displacement  there 
is  a  consequence,  rather  than  a  restriction. 

The  coefficients  and  are  evaluated  in  the  following  way.  A^^  is 

evaluated  by  considering  the  boundary  condition  at  x  =  0 

=  s  J  e  u^(o,t)  dt  (IIIj23) 

o 

The  coefficient  A^  must  be  formulated  along  the  characteristic  x  -  ct 

_s+ 

A2(s)  =  s  /  e  u^(ct,c)  dt  (lll.2lf) 

o 
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The  displacement  then  is 


u(x,t)  =  ^  (o, s)e^^"^^^^*ds  +  — ^  /  ^  (— , s)e^^”^2^ ®^*ds  (JH:*25) 


27ri  Br  27ri  Br  s 

If  this  expression  is  now  evaluated  for  large  s 


f^Cs)  =  - 


"  +  0 
2  ^2 


N_/ 

c  2dc^  d^s 

Thus  the  solution  along  x  +  ct  is 


u(x,t)  =  u  (t-  — —  t)  exp  ( — 

®  c  ^  2dc2 


(III. 26) 


x) 


1  +  0  (-  (t  -  -) 
a  c 


(111.27) 


Since  v  is  negative  with  respect  to  the  wave  front  this  solution  is  the  elastic 
wave  at  some  distance  from  the  end  of  the  rod.  The  attenuation  factor 

k  - 
^  ^  2 
2dc 

gives  the  amount  of  damping  in  the  elastic  wave. 

Inserting  representative  values  of  the  parameters  for  stress  loads  near  the 
yield  limit  for  the  metal  aluminum 


X  = 


a 


a  =  10  / 


cm 


p  =  2.7  gm/cra 


pJ)‘^ 

o  2 

D  =  10'^  cm  / sec 
V  =  oc  D  =  10^  cm/  sec 

-6 

the  average  distance  between  dislocations  d  is  taken  as  10  cm  and  the  sonic 
velocity  is  5-3  x  10^  cm/sec. 

Then  the  attenuation  coefficient  is 

10^  (10^) 


k  = 


(2)(2.7)(1^^)(5.3  X  105)^(10  +  6) 


=  6.5  X  10"  / 


cm 


(111,28) 
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which  is  in  reasonable  agreement  with  experimental  values  of  internal  friction. 

On  the  other  hand,  if  the  expression  for  u  is  integrated  for  large  t  away 
from  the  elastic  wave  front,  another  wave  front  and  the  boundary  layer  will  come 
into  prominence.  To  evaluate  the  wave,  consider  the  second  term  in  the  integral 
expression  for  u. 


(5) 


/  .\  1  r  fC  \  st-f_(s)x 

u(x,t)  =  ■; -  J  u  (— »s)  e  2'  '  I 


2Tt±  Br 


e  3 


ds 


(1X1.29) 


Rewriting  the  variable  in  the  exponential  as 


t  (s  -  A - ) 


u= -  J  u  (■?,s)e 

2'nl  Br  ®  " 


a. 


(III. 30) 


This  integral  may  be  evaluated  by  the  saddle  point  method  for  t  large.  The 
saddle  point  is  at  where 


1  ^2  (»2)  f  -  « 


The  prime  denotes  differentiation  with  respect  to  s.  The  dominant  term  in  the 
expansion  is 


u(x,t) 


(III.31) 


s^  2?r(x-x*)  t^(e} 


where  x*,  t*  are  the  x,  t  i)olnts  on  the  elastic  wave  front  corresponding  to 
this  wave, 

The  exponential  term  is  greatest  for  fixed  t  if  ~  From  the 

expression  for  this  occurs  at  Sg  =  0 
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Thus 


u(x,t)  = 


llm 

s->o 


u 


-  -  x*,t* 
c 


fg  (0)  =  2 


(c+v)(c-v)d 


3-1 
V  A, 


(111,32) 


u(x,t)  J - i  u  (ct,t)  dt 

V  47r(c  -V  )(x-x*)  O 

which  is  the  value  at  the  saddle  point. 

The  propagation  of  stress  waves  down  the  elastic  rod  is  seen  to  consist  of 
a  damped  elastic  wave  which  propagates  at  the  sonic  velocity  and  a  number  of 
dislocations  which  move  backwards  from  the  elastic  wave  front,  carryiiig  with 
them  additional  deformation.  The  boundary  conditions  at  the  end  of  the  rod 
which  are  initially  given  as  a  single  stress,  decompose  into  an  elastic  displace¬ 
ment  given  by  the  linear  stress- strain  law  and  an  added  displacement,  given 
mathematically  by  the  boundary  layer,  which  represents  permanent  plastic  strain. 
Nowhere  in  this  solution  is  there  a  forward  propagating  wave  that  moves  at  a 
velocity  lower  than  the  elastic  sonic  velocity. 


IV.  Plastic  Deformation  of  a  Semi- Infinite  Half  Space  Under  a  Suddenly  Applied 
Load  on  the  Surface 

In  this  section,  the  solution  of  the  plastic  and  elastic  deformation  of  an 
elastic  half  space  with  a  continuous  distribution  of  dislocations  will  be  studied. 
The  problem  will  he  restricted  for  mathematical  convenience  to  the  plane  strain 
problem.  This  corresponds  to  a  load  whose  area  of  application  extends  to  infinity 
in  the  one  direction  and  is  a  prescribed  function  of  position  in  another. 

The  equations  of  motion  are 

Sv, 


-11+  p-1 


m 


St. 


(IV.I) 


^11  =  ^  =11 


(IV.2) 


Da  d.  .  =  ^  d.  . 

Assuming  Hooke ' s  law  for  the  elastic  part 


(IV. 3) 


c.  .  =  i(u..  .  +  u. .  . ) 

ij  2'  i,j 


d,  •=  c.,  *"  c,,  6-. 

Ij  ij  kk  ij 


(iv.Jt) 

(IV.5) 

(IV.6) 


where  are  Lame's  constants 

Since  the  medlijm  undergoing  deformation  is  isotropic 

s  =  S6  .  +  t . 

im  im  Im 


(IV.7) 
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Restricting  the  equations  to  the  case  of  plane  strain,  the  equations  of 


motion  then  are 


V^u  -  (2-n-l) 


3x 


1  ^ 
fi  bx 


1  XX 
a  hx 


^  St 
1  xz 


,  M 


Sz 


i!!u 


V^w  -  (2r)-l)  ^  -  i  ^ 

Sx  /i  Sz 

.  St  ^  St  1  S^ 

1  xz  1  zz  1  w 


/Li.  Sx  p,  Sz 


5^ 


^11  =  I;  =11 


Da  d.  .  =  ^  d.  . 

St 


(IV.8) 


(IV.9) 


(iv.io) 

(IV.  11) 


Where  u,  w  are  the  displacements  in  the  x,  z  directions  respectively, 

■+■  =  n 

and  A  is  the  dilatation 

Su  ^  ^ 

Sx  Sz 

These  equations  take  on  a  more  instructive  form  through  the  application  of 
potential  theory. 

Introduce  the  potential  functions  0,  0,  S,  T  defined  by 


Civ.13) 


=  ^+  ^+  2-^+2^ 
Sz  dx  ^2 


IHfferentiating  the  first  equation  with  respect  to  z  and  the  second  with 
respect  to  x,  it  is  readily  observed  that  the  plastic  driving  stresses  are 
harmonic  functions 


t 


t 


t 


t 


XX 


yy 


xy 


M  - 

T)V^  0 

Sz 

^  _ 

dz 

ay^ 

^  - 

'd^f 

dz 

ax^ 

^y 

3x3y 

(IV.14) 


(IV.15) 


(IV.16) 

(IV.17) 


Inserting  these  values  into  the  equations,  the  equations  reduce  to  the 
following 

^  =  V^S  (IV.18) 


1 

2 

=  V^T 

(IV.19) 

(c^V^  - 

^  )  7^0  + 

(DgV^  ■ 

=  0 

(IV.20) 

t\j 

at 

(c^V^  - 

^  )  7^0  + 

(DxV^  • 

7^0  =  0 

(IV.21) 

'  i 

at 

where  c  ^ 
1 


P  2  P 


From  the  nature  of  these  equations,  four  houndary  conditions  must  he  pres¬ 
cribed  for  the  free  surface. 

The  first  two  equations  Indicate  the  propagation  of  two  elastic  waves.  The 
second  pair  of  equations  are  the  same  form  as  the  one  dimensional  plastic-elastic 
wave  studied  in  a  previous  section.  We  recall  a  theorem  derived  therelnj  the 
derivative  of  a  term  involving  the  dislocation  parameters  may  be  replaced  by 
division  by  the  average  spacing  of  the  dislocations 


D 

D  div  =  =  V 

1  k 


(IV.  22) 


where  V  is  the  dislocation  velocity.  Then  the  second  pair  of  equations  is  identical 
in  form  to  the  one  dimensional  case 


(c^V^ 
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)  ■¥ 


(V2  grad 


(IV.23) 


(c^V^ 


— —  )  +  (vi  grad  -  — )  V  0  »  0 

at^  at 


(IV. 24) 


The  nature  of  these  solutions  is  recapitulated  here.  One  part  of  the  stress 
wave  propagates  at  sonic  velocity,  the  wave  front  being  damped  by  the  action  of 
the  dislocations.  The  dislocations  move  away  from  the  wave  front  and  further 
displacement  occurs  behind  the  wave  front  and  at  the  same  time,  the  stress  state 
is  relieved. 

Here  the  only  difference  is  that  the  non-elastic  behavior  occurs  in  both  the 
shear  and  dilatational  modes.  Since  the  diffusion  coefficient  that  appears  in 
these  equations  is  dependent  on  the  maximum  shear  stress,  more  plastic  deformation 
is  to  be  expected  in  the  shear  wave. 
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A  detailed  solution  is  evaluated  for  a  simple  dependence  of  the  plastic 
driving  force.  In  this  simple  case,  the  tensor  components  are  all  assiimed  to 
vanish  with  the  exception  of  the  scalar  term. 

Under  this  assumption,  the  equations  of  motion  for  elastic  and  plastic  defor¬ 
mation  of  the  half  space  are 


V^u  -  (2r)^l)  ^  -  i 

Bx  /J.  Sx 


1_  ^ 


(IV.25) 


-  l2r,-l)  M  .  1  ^ 
Sz  ^  Sz 


1_  ^ 


ot 


(IV.26) 


(IV.  27) 


The  boundary  conditions  are 

^  =  P  (x,t)  z  =  0  t  >  0 

z  o  ~ 

0=0  z  =  0  ^>0 

T=0  z  =  0  t>0 

xz  — 

These  boundary  conditions  state  that  a  surface  traction  normal  to  the 
surface,  given  by  P^,  is  applied  to  the  surface,  that  the  surface  is  free  of  shear 
tractions,  and  that  the  driving  force  for  plastic  deformation  and  dislocation 
movement  vanishes  at  the  free  surface. 

These  equations  are  similar  to  the  equations  of  motion  for  deformation  of  a 
porous  media.  The  following  method  of  solution  of  these  equations  is  by  potential 
functions  and  the  transform  calcifLus.  The  method  of  solution  follows  very  closely 

(15) 

that  of  McNamee  and  Gibson  .  Since  the  thermodynamic  treatment  of  dislocation 
movement  indicates  that  the  dislocations  move  at  nearly  a  constant  velocity  when 
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the  shear  stress  is  above  the  static  yield  stress,  and  that  the  velocity  does  not 
depend  on  the  magnitude  of  the  stress  above  the  yield  stress,  little  dependence 
on  the  transient  inertial  forces  is  expected.  Therefore,  the  inertial  forces  are 
taken  as  zero  in  the  following  development  of  the  equation. 

^  =  0 

The  solution  of  these  equations  means  that  a  static  elastic 
established  instantaneously  after  the  load  is  applied.  Then  the 
matlon  ensues. 

The  equations of  motion  reduce  to 

V^u  -  (2r)^l)  _  1  ^  0 

hx  IJ.  dx 

V^w  *.  (27)-1)  M  _  i  ^  =  0 

hz  M 

D(t)V^(A)  = 

at 

%  differentiation  of  the  first  two  equations, 

(0+2^7^^)  =  0 

so  that,  if  S  is  a  harmonic  function  of  x  and  z 

a  =  2^  ,t]A)  (IV,33) 

az 

Another  potential  function  E  may  be  introduced  such  that 

A  =  V^E  (IV.  3^^) 


(IV.28) 

stress  field  is 
plastic  defor- 

(IV.29) 

(IV, 30) 

(IV.31) 

(IV.32) 
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The  displacements  u  and  w  are  related  to  these  potential  functions 


u  =  - 


3x 


(IV.35) 


(IV. 36) 


dz  dz 

Inserting  these  expressions  in  the  equations  of  motion,  the  potential 
functions  satisfy 


D  (t)  E  =V^ 


(IV. 37) 


S  =  0 


(IV. 38) 


Utilizing  the  linear  elastic  stress-strain  relations,  the  components  of 
stress  my  be  written  in  terms  of  the  potentials 


ox  _ 

-z 

3S 

2m  ' 

dz^ 

3x2 

3z 

az  _ 

-z  ^ 

3s 

2m  ” 

’  ax^ 

dz2 

3z 

^xz  _ 

d^E. 

z 

2m 

3xdz 

3x3  z 

(IV. 39) 


(lv.4o) 


(iv.4i) 


^  -  TjV^E  (IV  42) 

2/J.  dz 


These  expressions  identically  satisfy  the  equation  of  motion. 
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Since  the  coefficient  D  is  a  function  of  the  maximum  shearing  stress,  it  is 
more  convenient  to. solve  these  equations  in  a  non-dimensional  form.  The  substi¬ 
tutions  are 

X  =  bx' 
z  =  bz' 
t  = 

where  b  is  a  characteristic  length  dependent  on  the  original  load.  This  length 
will  be  clarified  later. 

This  substitution  introduces  a  non-uniform  time  scale  into  the  problem,  but 
will  not  introduce  any  serious  difficrolties. 

The  equations  for  the  potentials  become 

V^E«  =  /  III  (IV.43) 

2 

V  s‘  =  0  (IV. 44) 

The  superscript  primes  will  be  dropped  as  no  confusion  will  be  encountered.  The 

boundary  conditions  become 

a  =  P  (x,t) 
z  o  ' 

0  =  0  z  =  0t>0 

T  =0 
xz 

The  formal  solutions  of  these  equations  will  now  be  developed  by  the  integral 
transform  techniqSie. 

Introduce  the  Pourier-Laplace  transform  of  the  potentials 

00  00 

dt  /  cos  (^x)  E  (x,z,t)  dx  (IV. 45) 
o  o 
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(lV.lt6) 


S  (^>z,p)  =  /  e”^  dt  /  C08  €x)  S  (x,z,t)  dx 
o  o 

The  differential  rdatlons  of  this  transfom  are 

—  =  -p  E  (IV.47) 

dt 


(iv.W) 

(IV.49) 


vfaere  is  the  corresponding  Fourier  sine  transform. 


The  transfom  stresses  and  displacements  nay  he  written 
2/i  dz^  dz 


(IV.  50) 
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(IV.  51) 


(IV.52) 


(IV.53) 


(IV.  5^) 


(IV.  55) 
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The  differential  equation  for  the  transform  potentials  become 


(  IL  -  -P  )  (  ^  ^  ^  Q 


Sz‘ 


(IV,  56) 


,  s=o 

subject  to  the  boundary  conditions 
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T  =0 
xz 


u  -*n 


z  =  0 


t  >  0 


as  z  ->  ea 


(IV.57) 


The  .":oiier'0.  solutions  for  the  transform  potentials  may  be  written 

■  z(^^  + 


E  =  e 


(IV, 56) 


(IV.  59) 


The  coefficients  may  be  calculated  by  considering  the  values  of  the  stresses 
at  the  boundaries. 


For  a  =P^atz=0 
z  o 
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(IV. 60) 


For 


a  =  0  at  z  =  0 


A2'np  +  Bi^  =  0 


(IV, 6i) 


For 


T  =  0  at  z 
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.li  +  As  (e^+p)  =  0 


(IV.  69) 
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Thus 


A,  -  .  ^  _Slif±piV2 

(IV.  63) 

-Q 

(IV.64) 

/s 

r,  _  T)^p 

-  -  -sjr 

(IV.  65) 

Q  =  +  ?*  + 

(IV.66) 

the  change  of  scale 

P  =  4^8 

2/1 

(IV.67) 

A.,  "o 

2/i|2g 

(IV,68) 

8i  =  -  ^ 

2/i^G 

(IV.69) 

G  =  (l+T/s)  -  (1  +  s)^/^ 

(IV.  70) 

The  tilde  on  the  pressure  transform  indicates  that  the  transform  is  exposed 
as  a  function  of  the  variables  ^  and  s . 

The  potential  transforms  may  he  written 
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The  original  potentials  are  given  by  inverting  the  transforms 


E  = 


_  a 


,  / 

47r/il  Br 


ds  / 


j/  2 


COS  (^x)  d4 


(IV.  73) 


S  = 
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47r^tl 


/  e®^  ^ds  /  ■  °^®® -  cos  (4x)  d^ 

Br  o  G 


(IV. 74) 


where  Br  denotes  the  Bronwich  contour  in  the  coii5)lex  s-plane. 

The  components  of  stress  and  displacement  may  be  formulated  by  inserting 
transform  potentials  in  the  expressions  for  the  stresses  and  displacements,  per¬ 
forming  the  indicated  differentiations  and  inverting.  After  some  manipulation, 
these  are 


47T/il  Br 


/  e®^^%sr‘^o 


. 

G 


cos  ^x  d| 
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2/i  4ttp1  Br  o 
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cos  (|x)  <14 


For  later  convenience  in  discussing  the  dynamic  stress-strain  curves,  the 
formula  for  the  linear  strain dw/dz  is  included 
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These  Integrals  will  be  evaluated  for  a  particular  load.  The  load  is  a 

constant  pressiare  on  a  strip  of  width  2b  in  the  x  direction  on  the  surface  t  =  0. 

The  load  extends  to  infinity  in  the  y  direction.  This  load  is  ai^lled  suddenly 

on  the  surface  at  time  t  =  0.  Thus  the  pressure  is  represented  by  the  function 

(j=P*l(t)  x<b 

z  o 

0=0  X  >  b 

z 

where  1  (t)  denotes  Heaviside's  unit  function 
0  t  <  0 

1  (t)  = 

1  t  >  0 

If  the  dimensionless  coordinates  are  Introduced  with  b  as  the  characteristic 
length  the  pressure  is 


cw  w 

a„  (4^0,p)  =  /  e“^  dt  /  a  (x,0,t)  cos  (4x)  d|  (IV.S2) 


=  2  ^o 
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n  |®s 

Due  to  the  con5)lexity  of  the  ensuing  Integrals,  for  arbitrary  position,  they 

will  be  evaluated  here  for  z  =  0,  x  =  0.  By  the  boundary  conditions  o  =  P 

Z  Oy 

T  =  0.  By  invoking  symmetry  arguments  on  the  horizontal  displacement  u,  u  =  0. 
3CZ 


The  other  integrals  to  be  evaluated  are 
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with  t  = 


Since  these  integrals  are  convergent,  either  absolutely  or  in  the  Cesaro 
sense,  the  order  of  integration  may  he  changed.  The  integrals  may  be  written 
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where 
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The  evaluation  of  the  inverse  Laplace  transforms  will'  be  carried  ont  for 
the  explicit  elastic  parameter  ratio  rj  =  1.  This  corresponds  to  an  elastically 


incompressible  medium  with  Poisson's  ratio  v=  0.5 
The  integral  for  is 


rr 


Tri  Br  (l+s)  -  (l+s) 


ds 


This  complex  Integral  has  two  poles  at  s  =  -1  and  s  =  0  and  a  branch  point 
at  s  =  -1.  Rewriting  the  integral 
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Applying  the  Laplace  shift  theorem 
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The  integrals  for  the  vertical  displacement  and  horizontal  stress  become 
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These  integrals  with  respect  to  g  do  not  converge  absolutely.  However,  these 
integrals  converge  in  the  Cesaro  sense  and  represent  transform  images  of  fxmctions 

which  are  finite  for  finite,  t,  z  and  x.  To  demonstrate  the  convergence  in  the 

2Po  3^^ 

Cesaro  sense,  a  term  of  the  form - -  is  subtracted  from  each  integral. 

^r^l  g^ 

These  integrals  then  tend  to  a  finite  value.  These  integrals  then  exist  as  Cauchy 
principle  values  and  may  be  evaluated. 

1  ^  sJ^  (gt^^^)  dg  =  (7rt)^^%rf  — ^2-  Ei  (-  -i-  )  (iv.iCl) 
TT  o  g^  2t  ^  '41 

i  /  Mill  dg  =  (•n-t)^'^  erf  ( — \rT^  +  TT^  -Si  (~)  (iv.l02) 

TT  o  g^  2t  /  4t 
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where 


-  Ei  (-x)  =  ;  I  *  dt 


Si  (  x)  =  / 


sin  t 


dt 


Inserting  these  expressions  into  the  equations  for  the  displacement  and  stress 

(IV. 103) 
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^  =  erf  (-4/?)+  tV^  Si  (^) 

2t  /  4t 

^  ^  =  erf  (-^Ti)  +  Si  (-  — ) 

az  2tV^  2t 


(IV.104) 

(IV. 105) 


These  expressions  are  plotted  in  Figure  8.  The  maximum  shearing  stress,  by 
symmetry  given  as 

1 CT  a  I 

T  X  -  Z 

ma*  =  - ^ - >  (IV.106) 

is  also  plotted. 

The  time  scale  is  the  dimensionless  time  scale.  The  actual  time  scale  for 
alum1  n\im  is  plotted  under  the  following  restrictions,  above 

the  yield  shear  stress  =  10^  cm  /sec. 

For  short  times  after  the  impulse  little  plastic  deformation  has  taken  place. 
Around  100  jx  sec  the  plastic  deformation  becomes  appreciable  and  the  decay  of  the 
shear  stress  becomes  rapid.  The  horizontal  stress  rises  rapidly^  approaching  the 
value  of  the  vertical  stress.  If  the  shear  stress  is  allowed  to  decay  all  the 
way  to  zero,  then  the  horizontal  and  vertical  stresses  are  equal  and  a  state  of 
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i'^'lgure  8  -  SbrcBses  and  Strains  at  x  =  0  2  =  0 

Strip  Load  of  Width  2B  on  Ssai- Infinite  Halfspace 


I 
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Dirsensionless  Tine 


l^drostatic  stress  exists. 

The  stress- strain  curve  that  is  observed  during  deformation  under  rapid 
loading  is  not  a  \uii’que  curve.  Its  basic  conqjosition  consists  of  two  mechanisms, 
the  first  elastic  and  the  second  dislocation  motion.  The  first  mechanism  responds 
instantaneously  to  the  inqnilse  while  the  second  responds  at  a  given  rate  charac¬ 
teristic  of  the  medium. 

The  nature  of  this  stress- strain  relationship  may  be  reasonably  expected  to 
vary  at  different  points  within  the  medium,  since  the  original  stress  distribution 
will  vary,  and  thus  the  average  motion  of  the  dislocation  will  vary.  The  relax¬ 
ation  of  stress  within  the  medltim  is  of  course  dependent  on  the  movement  of  the 
dislocations. 

The  stress  and  strain  are  functions  of  the  dimensionless  time  variable 


a 


To  convert  from  the  reduced  time  to  actual  time,  the  dependence  of  D  on  the  max¬ 
imum  shear  stress  must  be  known.  This  was  derived  in  the  previous  section. 


D  =  T  > 

o  o 

where  is  the  attftic  yield  stress,  k  is  Boltzmann's  constant  and  T  is  the 

absolute  ten5)erature.  A  representative  diffusion  coefficient  is  shown  in  Fig.  9. 

Using  the  graphical  solutions,  the  observed  stress- strain  curves  were  numer¬ 
ically  calcialated.  These  are  shown  in  Figure  10  for  various  times  after  in^pact. 
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Diffusion  Coefficient  For  Dislocation  Motion 
Versus 

Maximum  Shear  Stress 


The  stress-strain  curve  Initially  is  a  straight  line  corresponding  to  pure 
elastic  deformation  at  very  early  times.  At  later  times  the  stress- strain  curve 
has  a  concave  and  a  convex  part.  This  is  explained  that  below  the  yield  stress, 
the  original  elastic  stress  strain  cuirve  is  followed.  Above  the  static  yield 
stress,  the  strain  component  increases,  but  to  a  finite  limit,  while  the  stress 
is  maintained.  At  points  away  from  the  boundary,  the  stress  will  change  as  a 
function  of  time  in  addition  to  the  strain  change.  It  is  concluded  that  the 
observed  stress-strain  curve  is  very  much  dependent  on  the  geometry  of  the  system 
and  where  the  relationship  is  measured. 

Examine  now  the  point  on  these  stress- strain  curves  where  the  first  con¬ 
cavity  appears.  A  sharp  turn  is  hard  to  describe  and  for  this  particular  set 
of  curves  no  actual  sharp  yield  point  is  indicated.  The  two  per  cent  offset 
cxirve,  as  good  a  definition  of  the  dynamic  yield  point  as  any,  is  plotted  as  a 
function  of  time  in  Figure  11. 
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V,  Impact  Plastic  Deformation  -  Cylindrical  Shell  Approximation 


The  equations  of  motion  for  dynamic  plastic  deformation  are  very  cumbersome 
when  applied  to  a  cylindrical  shelJ .  The  difficulty  is  caused  by  the  nature  of 
the  transcendental  integrands  that  occur  when  six  boundary  conditions  must  be 
satisfied. 

To  overcome  the  computational  difficulties,  a  shell  approximation  is  derived 
and  the  solution  to  the  resulting  set  of  equations  is  obtained.  The  establishment 
of  the  shell  equations  is  guided  by  the  elementary  three  dimensional  equations  pre¬ 
viously  derived ^ 

The  three  dimensional  solution  indicates  that  the  physical  parameters  in  the 
equations  of  motion  are  the  displacement  gradients^C which  may  be  resolved  into 
elastic  strains  and  permanent  displacement  gradient;.;  '>r  nJ.astic  strain)  the 
stresses  associated  with  the  elastic  strain,  and  a  stress  associated  with  the 
instability  of  the  elastic  strain  configuration, 

We  therefore  seek  a  formulation  of  the  problem  which  will  include  the  effect 
of  these  parameters.  To  accomplish  this,  a  variational  teclinlque  is  used. 

The  formulation  of  a  variational  principle  needs  expressions  for  the  kinetic 
and  potential  energy.  No  rigorous  proof  shall  be  given  for  the  existence  of  such 
potentials.  Their  justification  is  based  on  an  appeal  by  analogy  to  the  results 
obtained  in  the  three-dimerif’.ion.u  I  c-i.-.e. 

As  in  the  previous  solutions,  the  inertia  terms  will  be  omitted.  The  stress 
wave  propagation  across  the  thickness  of  the  shell  will  be  assumed  to  have  occurred 
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many  times  and  a  quasi-steady  state  situation  to  be  established. 


The  potential  energy  of  the  deformed  shell  consists  of  the  elastic  strain 
energy 


2 


2  ^ij 


(V.l) 


where  is  the  total  stress  and  is  the  geometric  strain  plus  the  potential 
energy  due  to  the  plastic  driving  force  minus  that  plastic  energy  which  is 
dissipated  by  the  diffusive  flow  of  dislocations.  Since  the  plastic  driving  force 
in  the  three-dimensional  equations  is  satisfied  by  a  potential  function,  it  is 
assumed  that  the  potential  energy  associated  with  it  may  be  written  analogously 
to  the  elastic  strain  energy.  The  gradient  of  the  plastic  energy  dissipates  itself 
by  the  diffusive  mechanism.  Mathematically  this  may  be  formulated 


2  n  t 

V  =  2-i-£ - -  t  (v  )dt 

^  •  2M  2M  ■' 

o  * 


(Vv2)' 


where  M  is  a  stress  modulus  which  is  taken  as  the  yield  stress  in  compression  and 
D  is  the  shear  stress  dependent  diffusion  coefficient  for  movement  of  the  disloc¬ 
ations.  A  simple  behavior  of  D  with  the  maximum  shear  stress  will  be  assumed.  If 
the  shear  stress  is  greater  than  the  static  yield  point 


D  =  D 

o 

If  the  shear  stress  is  lower  than  the  static  yield  point 


D  =  0 


Consider  an  infinite  cylindrical  shell  of  thickness  h  whose  middle  surface 
radius  is  R  (Figure  12),  The  coordinate  system  is  cylindrical  polar.  The  form¬ 
ulation  of  the  present  equations  is  restricted  to  the  case  where  there  is  no 
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dlsplacanent  or  variation  of  stresses  in  the  direction  of  the  generators,  i.e., 

u  =  o^=o.  (;v.3) 

Using  the  classic  approach  to  shell  theory,  the  only  strain  term  that  is  considered 
is  the  strain  along  the  middle  surface  in  the  tangential  direction  The  strain 

in  the  radial  direction  is  assvnned  to  be  negligibly  small,  as  is  the  shear  strain. 


The  linear  approximation  to  the  elastic  strains  and  stresses  will  overestimate 
the  elastic  displacement.  This  error  is  not  serious  in  the  computation  of  the 
plastic  strain  as  the  main  influence  of  the  elastic  stress  state  is  to  drive  the 
plastic  deformati.on  at  a  given  rate.  The  primary  information  needed  about  the 
elastic  stress  is  to  know  when  this  stress  exceeds  the  static  yield  condition.  The 
strain  is  expanded  as  a  function  of  the  strain  of  the  middle  surface  plus  a  linear 
term  in  r  (7) 


^1  ■  2  ^ 

=  . 


1  RUe 


0 


/-V 


(V.5) 

(V.6) 


^0" 

where  v  is  the  tangential  displacement  and  w  is  the  radial  displacement, 
r^resents  the  membrane  or  extensional  mode  of  deformation  while  X  ^r  represents 
the  strain  due  to  bending. 

The  elastic  stress  is  given  by 


(V,7) 


where  E  is  Young's  modulus 
and  V  is  Poisson's  ratio 
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The  total  stress  is  the  sum  of  the  elastic  stress  and  the  plastic  stress. 


a  = 


E 


^  2  0 
1-v 


P 


(V.8) 


The  plastic  driving  force  gradient  is  assumed  constant  across  the  thickness  and 
the  radial  gradient  then  is  zero. 


Pi,  i  -  R  ^Pe  P^ 


(v.p) 


The  potential  energy  may  then  be  written  using  subscripts  to  denote 
lifferentiation. 


V  = 


E 


r„  27T 

/  ^  / 


2(1-1^)  o 


(ve  ‘w)"  ^2 


R 


4r" 


Rdrde 


,27T 


-  /  2  / 


"l  ° 


^  f- 


(Vg  p(v^^wj 


Rdrde 


2  1 

“  (P0  +  P)^  ^  JRdrd© 


2R^M  o 


/•27r 

=  -  /  FwRd0 


(v.io) 


where  F  is  the  applied  pressure  in  the  inward  radial  direction. 


Carrying  the  indicated  integration  with  respect  to  r 


E  r2nl  h  ,  n2  h-^  ,  x2 

- p~  /  ~p  -  w)  +  - j-  (vq  +  w  J 

2(1-1^)  o  [R  ^  12E^  ^ 


Rd0 


-  (ve  -  w)  +  ^  P  (vg  +  Wg^ 

^  ,27r('  Dh  f  ^  x2  hp^ 

+  J  J  (P  +  Pa  )  dt  - 

o  !2R%  o  fc'  2M 

,2 

=  -  /  FwRde 


Rd  0 

Rd  0 


(V.ll) 
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The  differential  equations  of  motion  are  evaluated  by  minimizing  this  integral. 
Applying  the  principles  of  the  calculus  of  variations  this  is  straight  forwai^  and 
the  resulting  equations  are 


ft  h  ^  Dh^  Dh  ^  Dh^  ,  h  - 

R  ^ee-^  ''e  +  -3*'6Be-^Pe -MP0  =  o 

12R  R 


ft  h  ^  l)h^  ft  h  ,  Dh^ 

^  R  ^  Q  ^fYY)n  ^  R  ^ 

W  12rJ  t**7C7  l\  u  22Ir 


^  Pee  ^  ^  (^-^3) 


-  h  (vg  -  <r)  -  ^  (ire+  Wgg) 


*  ®  I  (ffle  *  *  f>  **  -  T  ■  ° 


('V.i4) 


with  D  = 


E 


(1-  v2) 


Introduce  now  the  dimensionless  displacements 


(V.15) 


w^  *  H 
R 


and  the  dimensionless  time 


(V.l6) 


The  equations  may  now  be  written.  The  third  equation  is  also  differentiated 
with  respect  to  time 

“i  <’'00  -  'e  '  *  ”2  <'''@0 ■"  ’''000>  -<”3  *  V  »'0  “  ° 

\  (w*  -  V'  0)  +  i^’eeee  +  ^'eee  )  +  ^3  p'  "  P'ee  +  p’  =  0  (v.i8) 

"®3  et  ■  '^'t^  ‘  \  et  ^  ^  ^P’ee-"  P't  =  ° 
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where 


°2  =  °1 

12R 


Eh 


12(1- v^)R^ 


D3  =  1 

”6-5 


The  primes  will  now  be  dropped  as  no  confusion  can  arise. 


Solution 


The  solution  to  these  equations  will  be  sought  in  the  transform  space. 
Denote  the  laPlace  transform  of  the  variables. 


w  (s,  0)  =  / 

e  w  (t, 0 )  dt 

(V.20) 

0 

^  (s,  =Tq 

e  V  ( t,  ®  )  dt 

(V,21) 

p  (s,  =  r 

e  p  (t,  0)  dt 

(V.22) 

o 


Assume  that  the  transform  variables  are  sparable  and  e]q>and  them  in  a 
Fourier  series. 

00 

w  (s,  0)  =  2  w  (s)  cos  nd  (V.23) 

^  n=0  n 

00 

V*  (s,  0)  =  v^  (s)  sin  n0  (■V.24) 

(s,  0)  =  p^  (s)  cos  n®  (V.25) 
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utilizing  the  transform  relations  for  derivatives,  the  equations  in  the 


transformed  variables  may  be  written 

(D^  +  D^)  +  (dD^  +  n^D^)  ° 

(-nD,  -  n^D^)  v  +  (D  +  n**^D„)  w  +  (D-  +  n^D,,)  p  =  f 
'1  3  n  1  2'  n  '  3  4'  ^n  n 


(V.26) 

(V.27) 


-ns  (D^  +  Dj^)  -  s  (D^  +  n^Dj^)  w^  +  (sD^  -  (n  +  l)^D^)  =  0  (V.28) 


where 

OO 

F  (9,0)  =  Z  ^  (s)  cos  ne  (V.29) 

is  the  liiPlace  Fourier  transform  of  the  applied  pressiore.  The  three  equations 
are  now  reduced  to  a  set  of  algebraic  equations  in  the  transform  Fourier  coefficients. 
Using  Cramer's  rule 

I,  ^  „2  +  D,.^2  +  _2  D,,  +  n  (nfl)^  D^)] 

.30) 


Vrx  cuucx  o  X  axe 

[  2  /D^  +  D-  D,-»  ^  2  .D^  +  D,  x2  +  2  D,.  +  n  (nH)^  D.  D,.)] 

(?  =  -  f  L-n  (  1  6  2  3)  s  +  n  (  3  ^^)  n  (  1  ^  2  - 

nj.  g  ^ 

** 


yK 

V 


^  (n  Dg  +  n^D2Dg)s  -  n^  +  n^  (nH)"^ 

-5  Y  A  “  ~ 


(V.31) 


where 


A  = 


A  = 


For  shells 


-n^  +  D^)  nD^  +  n-'^Dg 

3 

-nD^  -  n  Dg  +  n 


n(D^  +  Dj^) 


D3  +  n^Dj^ 


-ns(D3  +  Dj^)  -  s(D3  +  n^Dj^)  (sD^  -  (n  +  1)^0^) 
[(n^+l)D^D^^  +  (n^-l)D^D2Dg  +  "  (n^+l)  02^3" 

-2n^D2D3D,^]  s  -(n+l)^  (n^-l)^ 


Di>D3 


and  for  most  materials 


(V.32) 


(V.33) 


D2»D3>d^ 
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Thus  the  e]q>resslons  for  the  displacements  may  be  simplified 

.  ”Vj 

”  (n^l)  D^DjDg 


-f 


n(Dj^Dg)s  - 


“  (n  -DD^D^Dg  (s-7^) 


where  7  =  (h+l)^ 

n 


(V.34) 

(V.35) 

(V.36) 


These  are  the  general  transform  solutions  for  the  Fourier  coefficients  for  an 
arbitrary  applied  pressure.  A  special  case  will  now  be  examined. 

Let  us  assume  that  a  radial  pressure  of  uniform  magnitude  is  applied  over 
a  region  of  the  cylinder  surface  and  zero  is  outside  that  region  (Figure  lj{)* 


F  =  F(t)  -  a  ^  0  ^  a 
=  0  elsewhere 

The  Fourier  cosine  expansion  of  this  pressvire  is 


F(e,t)  =  F^(t) 


00  .  1 
a  .  2  ^  sin  na  ^  i 

—  +  —  2  -  cos  n  0 

n=l  j 


•(V,37) 


i(V.,38) 


The  time  dependence  of  the  pressure  on  the  surface  will  be  a  decaying 
exponential 


F  =  F  e 
o 


'  /3t 


The  LaPlace  transform  of  the  series  is  then  given  by 


Hs, 


a  +  ^  2 


sin  ntt 
n 


cos  n  Q\ 
j 


I  7T  TT  n=l 

The  transforms of  the  Fourier  coefficients  then  are  given  by 
n  >  1 

2F  S 

n  ^  +  o  (ns  -  n~^)  sin  nOi- 

(s-7n)(s-/3)(n2-l) 


(V.39) 


i(W.4o.) 


( .V..4  ]) 
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2x10  '  3x10 


■  (iV.42) 


2¥  ,  3x  . 

P  ^  _  o  (s-n  )  sin  na _ 

^  (s-7j(s-/B)(n^-l) 


Inverting  these  integrals 

c+ioo 

V  {.t)*-  - — .  /  w  (s)  e^  ds 

n  '  2  TTi  .  n'  ^ 

C-19® 

c+i"  . 

V  (t)  =  .  /  V  (s)  ds 

n  '  27r  1  ,  n'  ^ 

c-loo 


(V,43) 

(V.U4-) 


hy  the  values  at  the  residues  of  the  integrals 


»2  ,  (t)  .  g-tf” 

F  -)T(n^-l) 

o 


\  V  (t)  =-2_3in_na_ 
°  7r(n  -l) 


t  -1 


^-7 


/3-7 


n 


2  2  »  2 

/fi  -n  x \  7  "7  \ 

(.^  fa)e  ^  (  'n  'n)  e  ^ 


^-7 


A-7, 


(V,45) 


(V.46) 


These  results  vill  he  regrouped  to  indicate  the  elastic  and  plastic 
components.  To  show  how  this  is  done,  consider  the  case  of  the  Heaviside  unit 
pressure  pulse  i.e.,/3=  0 


\  ^  (t)  =  ga  sin.na 

F^  “  7r(n^-l) 


2  ,  2  . 
-HTn  (“7„-7„) 
- + -  e  n 


(V.Jv?) 


or 


!  ^\  2n  sin  na  f  ^7„t  2  ,, 

_2  w  (t)  - - 5 -  -  e  n  -  n  (l-e  n  ) 

F^  “  (n^-1)  L  J 


(y.48) 


At  the  instant  of  applied  load  no  plastic  deformation  is  present.  Only 
elastic  deformation  can  account  for  the  initial  displacement.  At  later  times  the 
elastic  deformation  decays  and  the  plastic  deformation  increases  as  the  dislocat¬ 
ions  have  a  chance  to  move.  In  the  above  expression  the  first  term  has  its 
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maximum  value  at  t  =  0,  and  thereafter  decreases  while  the  second  term  starts 
at  zero  and  Increases  with  time.  Thus  we  identify  the  first  term  with  the 
elastic  displacement  and  the  second  term  with  the  permanent  or  plastic  displace¬ 
ment,  (Figure  l4) . 


Taking  a  finite  value  of  the  same  regrouping  may  be  made. 


D  / . \  2n  sin  na 
2  w  (t)  =  - - - 

T  ^  TTCn-^-l)  L 


^2  V  (t)  =  -  ^ 

/  “  ^  (n^-1) 


-y  e  y  n~.  -  ft  -7  i;v 

_ n  (e  -e  n  J 


^  -y 
n 


A  -X 


n 


jB  -7 


n 


(V.49) 

(V.50) 


The  solution  for  the  zeroth  mode  requires  special  handling.  This  mode 
corresponds  to  a  xmiform  pressure  applied  around  the  circumference.  The  transform 
equations  for  the  radial  displacement  and  the  plastic  driving  stress  for  this  mode 
are 


D 

1 


w 

o 


D  p 

3  0 


-sD  w  +  (sD  -D  )p 
3  0  6  5  0 


=  0 


The  solution  for  w  is 
o 


w  = 

o 


F  (sD  -D  ) 
o  6  g 

D  (sD  -D  )+D^ 
1  6  5  3 


Inserting  the  value  for  ^ 


(V.51) 

(V.52) 


(V.53) 


(V.54) 
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Inverting  and  using  the  fact  that 

5  o 


where 


D  w 


~  ( t)  =  i 


.  ,  -(l+  6)t 

p-1-5  ^-1-  i 


hi. 


6  = 


(V.55) 


For  the  particular  example  to  be  evaluated  this  term  is  small  compared  to  the 
hi^er  modes  and  will  he  neglected. 


The  first  mode  requires  special  handling  for  solution.  This  mode  corresponds 
to  the  translation  of  the  whole  shell.  No  deformation  is  involved  and  therefore 
this  mode  is  not  evaluated. 


We  will  now  investigate  the  qualitative  behavior  of  the  deformation  modes. 

The  elastic  portion, starts  at  some  initial  value  and  decays  to  zero  in  a  finite 
time.  The  plastic  component  starts  at  zero  displacement  and  increases  to  a 
maximum  at  the  same  time  the  elastic  part  reaches  zero.  The  solutions  beyond 
this  time  are  discarded  as  they  imply  that  the  plastic  or  permanent  deformation 
does  work  to  increase  the  elastic  deformation.  However  the  plastic  part  indicates 
an  amovint  of  energy  which  is  permanently  lost  to  the  recoverable  portion  of  the 
total  energy  of  deformation.  The  solutions  after  this  critical  time  indicate 
that  this  energy  is  recovered.  Since  we  have  made  the  identification  of  this 
displacement  with  nonrecoverable  deformation,  the  solution  beyond  the  critical 
time  is  meaningless. 

The  plastic  deformation  will  increase  imtil  the  maximum  shear  stress 
decreases  to  the  static  yield  point.  At  this  time  the  plastic  deformation  stops 
at  that  value.  This  may  be  shown  by  computing  the  actual  time  scale  from  the 
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h  =  2  cm 

R  =  100  cm 

M  =  10  Kbars 

E  =  750  Kbars 

o  p 

D  =  lO'^  cm  /  sec 

A  uniform  pressure  of  100  Kbars  is  applied  over  a  6o°  arc  on  the  exterior  of  the 
shelly  (Figure  iJ+a)  so  that  oc  =  30".  The  time  dependence  of  the  load  is  a  decaying 
exponential  (Figure  l4b) 

P  =  F  e”^^ 
o 

In  the  actual  time  scale  J3  is  chosen  as  10^/sec.  In  the  reduced  time  scale 
jS  =  10^ 

The  elastic  displacement  from  the  zeroth  mode  is  negligible  compared  with 
the  displacements  of  the  second  and  third  modes. 

The  first  mode  does  not  involve  any  distortion  of  the  shell  and  is  therefore 
neglected. 

We  retain  only  the  first  two  modes.  The  two  modes  n  =  2,3  were  computed. 

The  elastic  portion  of  the  solution  for  these  two  modes  was  calculated  first 
(Figure  I5  and  16).  The  maximum  shearing  stress  was  calculated  at  a  number  of 
points  to  determine  the  cutoff  time,  i.e.,  the  time  at  which  the  diffusion 
coefficient  went  to  zero.  The  plastic  portion  of  the  displacement  for  these  two 
modes  was  calculated  (Figure  I7  and  I8),  and  evaluated  at  the  cutoff  time.  Table 
I  shows  the  values  of  the  cutoff  points,  the  values  of  the  plastic  deformation  at 
that  time,  and  the  computed  values  of  the  radial  and  tangential  displacements. 
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value  of  the  dimensionless  time 


t 


(v.56) 


When  the  shear  stress  is  below  the  static  yield  value  t  =  00. 

To  compute  the  value  of  the  plastic  strain,  the  time  dependence  of  the  maximimi 
shear  stress  must  be  known. 


Now 

(V.5T) 

’■■M  Hr" 'el 

(V.58) 

where  ^  is  the  shear  modulus.  In  the  shell  approximation,  the  radial  strain  is 
assiomed  zero  so  that 


U 


(V.59) 


For  very  thin  shells,  the  tangential  strain  is  given  by  the  membrane  strain  for 
low  modes 


-  w  (y.60) 

For  each  Foxtrier  component 

Cq  cos  n0  =  (n  cos  n  6  (V.6l) 


But 


nv 

n 


w 

n 


(V.62) 


Thus 


T 

2M 


Z  w  cos  n0 
n=o  “ 


(V.63) 


A  niomerical  example  for  the  permanent  defonnationof  a  cylindrical  shell 
under  a  high,  rapidly  decaying  load  is  given.  For  this  example,  a  shell  with  the 
following  geometrical  and  physical  properties  is  chosen 
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R  =  100  cm 


h  =  2  cm 

P^=  100  kbars  =  E/6 
|3  =  10^/  sec. 

Figure  19  shows  the  shape  of  the  deformed  shell. 

Since  this  computation  retained  only  two  of  the  terms  of  a  slowly  converging 
Fourier  series,  the  results  show  a  great  deal  of  wohhle,  particiilarly  on  the 
backside  of  the  plate.  Estimating  the  value  of  the  elastic  strain  to  30  modes, 
the  plastic  deformation  term  is  negligible  bgyond  6  =  45“.  Setting  all  of  these 
permanent  displacements  equal  to  zero  and  retaining  only  the  two  mode  approximation 
to  the  plastic  deformation  results  in  a  more  physically  realizable  picture.  This 
deformation  pattern  is  shown  in  Figure  20. 
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Figure  20 

Permanent  Ifeforniation  of  a  Cylindrical 
Shall  Under  Time- Dependent  Pi’essure  Pulse 
Refined  Approximation 


VI  Svmmiary  and  Conclusions 


The  prolDlem  of  plastic  deformation  under  time  dependent  loads  has  been 
investigated  from  the  viewpoint  that  at  least  two  mechanisms  work  simultaneously 
in  this  process  of  deformation.  In  addition  to  the  usual  elastic  deformation, 
further  microscopic  deformation  by  a  specific  microscopic  mechanism  has  beer 
studied.  This  mecha.nism  is  the  motion  of  dislocations  under  an  applied  stress. 

A  set  of  equations  of  motion  for  a  deforming  body  were  derived  \mder  these 
ass\miptions  of  this  physical  model.  These  equations,  when  restricted  to  hydrostatic 
variation  of  the  non-equilibrated  driving  force  are  Identi'-al  with  the  equations 
of  motion  for  a  compacting  porous  medium.  Formal  solution  of  the  equations  of 
motion  for  a  propagating  stress  wave  where  the  material  is  stressed  above  the 
elastic  limit  indicates  that  a  damped  elastic  wave  propagates  at  the  usual  sonic 
velocity  while  further  deformation  occurs  near  the  boundary  due  to  diffusive  flow 
of  the  dislocations. 

A  detailed  analytical  and  numerical  solution  for  two  problems  involving 
plastic  deformation  was  carried  out.  The  first  problem  solved  was  the  deformation 
of  a  semi-infinite  half  space  under  a  strip  load.  The  displacements,  stresses 
and  total  strains  were  analytically  determined  in  integral  form  and  were  evaluated 
for  a  specific  point.  Several  interesting  conclusions  were  demonstrated. 

1.  The  stress-strain  curve  tliat  is  observed  is  a  function  of  time  duration 
of  load,  level  and  rate  of  applied  stress  and  the  geometry  of  the  problem. 

2.  The  dynamic  yield  stress  is  explained  simply  as  the  first  noticable 
deviation  from  the  pure  elastic  curve.  The  dynamic  yield  stress  is  a  result  of 
the  process,  rather  than  an  explicit  material  property. 

9k 
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The  equations  of  motion  were  applied  to  deformation  of  a  cylindrical  shell. 
A  shell  approximation  including  the  effects  of  the  dislocation  movement  were 
established  and  a  niomerical  approximation  calculated. 
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